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Abstract 

A review of the possible relationship between 7-ray sources and supernova remnants (SNRs) is pre- 
sented. Particular emphasis is given to the analysis of the observational status of the problem of cosmic 
ray acceleration at SNR shock fronts. All positional coincidences between SNRs and unidentified 7-ray 
sources listed in the Third EGRET Catalog at low Galactic latitudes are discussed on a case by case 
basis. For several coincidences of particular interest, new CO(J=1-0) and radio continuum maps are 
shown, and the mass content of the SNR surroundings is determined. The contribution to the 7-ray flux 
observed that might come from cosmic ray particles (particularly nuclei) locally accelerated at the SNR 
shock fronts is evaluated. We discuss the prospects for future research in this field and remark on the 
possibilities for observations with forthcoming 7-ray instruments. 
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Figure 1: Plot of source locations (galactic coordinates) in the COS-B catalog. The undashed region of 
the sky was surveyed by the satellite. Most of the sources had fluxes higher than 1.3 x 10"^ photons cm^^ 
s~^ above 100 MeV. From Bignami & Hermsen (1983). 



1 Introduction 

Gamma-ray astronomy has unveiled some of the most exotic and energetic objects in the universe: from 
supermassive black-holes in distant radio galaxies to radio-quiet pulsars and the still enigmatic gamma-ray 
bursts. However, it has been conspicuously less successful in achieving one of its original goals of shedding 
light on the sources of Galactic cosmic ray nuclei. In this report we focus on the remnants of galactic 
supernovae, and their possible association with discrete sources of (> 70 MeV) 7-rays, as seen by the En- 
ergetic Gamma-ray Experiment Telescope (EGRET). In doing so, we attempt to lay a framework in which 
the long-standing question of the supernova remnant origin of Galactic cosmic rays may be addressed. 

The first firm detection of celestial high-energy 7-rays was achieved by Clark, Garmire and Kraushaar 
using the Orbiting Solar Observatory (OSO-3), when they discovered that the plane of the Galaxy was a 
source of photons with E > 70 MeV (Clark, Garmire &: Kraushaar 1968; Kraushaar et al. 1972). Higher 
spatial resolution studies made with the SAS-2 satellite, launched in 1972, revealed individual sources of 
7-rays from the Vela pulsar (Thompson et al. 1975), and confirmed the high-energy emission from the 
Crab (Kniffen et al. 1974). The long Ufe of ESA's COS-B satelhte (1975-1982) produced another major 
breakthrough in 7-ray astronomy: for the first time a significant number of 7-ray sources were seen which 
could not be identified with objects known at other wavelengths (see Bignami &: Hermsen 1983, for a 
review of COS-B results). Figure ^ shows the region surveyed by COS-B and the point sources discovered, 
as reported in the second COS-B Catalog (Hermsen 1981, Swanenburg et al. 1981). 

In 1991, the EGRET telescope was launched onboard the Compton Gamma-Ray Observatory (see 
Gehrels & Shrader 2001, for a recent review). The Compton satellite (1991-2000), the heaviest orbital sci- 
entific payload at the time of its launch, had three other experiments apart from EGRET. All of them have 
contributed to our understanding of the 7-ray sky, although we shall particularly focus on EGRET results 
in this report. The Third EGRET (3EG) Catalog, whose point-like detections are shown in Figure [21 is 
now the latest and most complete source of information on high-energy 7-ray sources. It contains 271 
detections with high significance, including 5 pulsars, 1 solar flare, 66 blazar identifications, 1 radio galaxy 
(Gen A), 1 normal galaxy (LMC), and almost two hundred unidentified sources, ^ 80 of them located at 
low galactic latitudes (see Grenier 2001 and Romero 2001 for recent reviews). 

The detection of pulsed high-energy emission from some 7-ray sources, on one hand, and the identi- 
fication of Geminga as a radio quiet pulsar, on the other, have prompted several authors to explore the 
possibility that all unidentified low-latitude sources contained in earlier versions of the EGRET Catalog 
(i.e. the Second EGRET -2EG- Catalog, Thompson et al. 1995, 1996) could be pulsars as well (excepting 
a small extragalactic and isotropic component which should be seen through the disc of the Galaxy). In 
particular, Kaaret & Cottam (1996) used OB associations as pulsar tracers, finding a significant positional 
correlation with 2EG unidentified sources. A similar study, including SNRs and HII regions (considered 
as tracers of star forming regions and, hence, of possible pulsar concentrations) has been carried out by 
Yadigaroglu & Romani (1997), who also concluded that the pulsar hypothesis for the unidentified 2EG 
sources was consistent with the available information. 
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Figure 2: Plot of source locations (galactic coordinates) in the Third EGRET Catalog. Different popula- 
tions are marked in different grades of light colours. The size of the dots gives a qualitative idea of the 
detected flux. Prom Hartman et al. (1999). 



However, spectral analysis by Merck et al. (1996) and Zhang & Cheng (1998) showed that several 2EG 
sources were at odds with the pulsar explanation; the spectra of many sources are too different from what 
is expected from outer or polar gap models of pulsar emission. Time variability in the 7-ray flux of many 
sources (discussed below) also argued against a unique population behind the unidentified galactic 7-ray 
sources. 

Most likely, the unidentified 7-ray sources at low galactic latitudes are related to several different galac- 
tic populations (e.g. Grenier 1995, 2000; Gehrels et al. 2000; Romero 2001). Among them there are surely 
several new 7-ray pulsars (e.g. Kaspi et al. 2000, Zhang, Zhang & Cheng 2000, Torres, Butt & Camilo 
2001, Camilo et al. 2001, D'Amico et al. 2001, Mirabal et al. 2000, Mirabal & Halpern 2001, Halpern et 
al. 2002). Pulsars remain as the only confirmed low-latitude population, since pulsed 7-ray radiation has 
been already detected for at least six different sources (Thompson et al. 1999, Thompson 2001). Other 
populations might include X-ray transients (Romero et al. 2001), persistent microquasars (Parcdes et al. 

2000, Grenier 2001, Kaufman-Bernado et al. 2002), massive stars with strong stellar winds (Benaglia et al. 

2001, Benaglia & Romero 2003), isolated and magnetized stellar-size black holes (Punsly 1998a,b; Punsly 
et al. 2000), and middle- mass black holes (Dermer 1997). Finally, there is a possibility that some 7-ray 
sources could be generated by supernova remnants (SNRs), especially those interacting with, or located 
close to, molecular clouds (e.g. Montmerle 1979; Dorfi 1991, 2000; Aharonian, Drury & Volk 1994; Naito 
& Takahara 1994, Combi & Romero 1995; Aharonian & Atoyan 1996; Sturner, Dermer & Mattox 1996; 
Esposito et al. 1996; Combi et al. 1998, 2001; Butt et al. 2001). This review is devoted to discuss this 
latter possibility in light of recent observations. 

SNRs are thought to be the main source of both cosmic ray (CR) ions and electrons with energies below 
the knee in the galactic CR spectrum, at ~ 10^^ eV - however, see Plaga (2002) for alternate theories. The 
particle acceleration mechanism in individual SNRs is usually assumed to be diffusive shock acceleration, 
which naturally leads to a power-law population of relativistic particles. In the standard version of this 
mechanism (e.g. Bell 1978), particles are scattered by magnetohydrodynamic waves repeatedly through 
the shock front. Electrons suffer synchrotron losses, producing the non-thermal emission from radio to 
X-rays usually seen in shell-type SNRs. The maximum energy achieved depends on the shock speed and 
age as well as on any competing loss processes. In young SNRs, electrons can easily reach energies in excess 
of 1 TeV, where they produce X-rays by synchrotron mechanism (see, for example, Reynolds 1996, 1998). 
Non-thermal X-ray emission associated with shock acceleration has been clearly observed in at least 11 
SNRs, and this number seems to be steadily increasing with time. In the case of the very nearby remnant 
RX J0852. 0-4622 (also known as Vela Jr.) the discovery was originally made at X-rays (Aschenbach 1998) 
and only then the source was detected at radio wavelengths (Combi, Romero & Benaglia 1999). 

As early as 1979, Montmerle suggested that SNRs within OB stellar associations, i.e. star forming 
regions with plenty of molecular gas, could generate observable 7-ray sources. Montmerle himself provided 
statistical evidence for a correlation between COS-B sources and OB associations. Pollock (1985) presented 
further analysis of some COS-B sources in the same vein. Statistical correlation studies of EGRET sources 
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Figure 3: Distribution of the Green's SNRs (circles) together with EGRET unidentified sources (diamonds), 
shown in galactic coordinates. Some of the coincident pairs that are studied in this report are marked. 
The top marks are SNRs for which TeV radiation has been detected. From Mori (2001). 



and SNRs have been presented by Sturner & Dermer (1995), Sturner et al. (1996), Yadigaroglu & Romani 
(1997), and Romero et al. (1999a). These studies show that there is a high-confidence correlation between 
remnants and 7-ray sources. Figure 01 shows the distribution of the SNRs in Green's Catalog (2000) along 
with the 3EG unidentified sources. Some of the coincident pairs that are studied in this report are marked. 
It is worth noticing, however, that the EGRET observations may contain important information about 
the hadronic component of cosmic rays concerning only the low-energy domain (typically less than several 
tens of GeV). Hence, these observations alone cannot solve the problem of galactic cosmic rays. The lat- 
ter effort require TeV observations, and that is why we also discuss future TeV observations of EGRET 
sources in this review. Note, in passing, that whereas an all-sky EGRET map is available, only pointed 
TeV observations are possible. 

SNRs can produce high-energy 7-rays through nucleus-nucleus interactions leading to vr^-production 
and subsequent 7-decays. The resulting 7-ray luminosity will depend on the local enhancement of the CR 
energy density as well as on the density of the ambient media. However, as shown below in Eq. (|10|) . the 
expected fluxes of vr'^-decay gamma-rays from a SNR are generally well below the EGRET sensitivity (quite 
importantly, the estimate given by Eq. (|10|) is almost independent of the proton spectrum). Thus, any 
detection of gamma-ray flux from a supernova remnant neighborhood would imply a significant enhance- 
ment of gamma-ray production. Such an amplification would be possible only through the assumption of 
the interaction of CRs accelerated by the shell of the SNR with a nearby high density environment - e.g. 
a dense molecular cloud. GeV 7-rays (and also TeV photons, see e.g. Pohl 1996) can be produced also 
by electrons through relativistic Bremsstrahlung and inverse Compton upscattering of cosmic microwave 
background photons, diffuse Galactic infrared/optical radiation, and/or the radiation field of the remnant 
itself (e.g. Mastichiadis 1996, de Jager & Mastichiadis 1997). Gaisser et al. (1998) modeled these processes 
in detail in order to fit the observational data for the SNRs IC 443 and 7-Cygni. Sturner et al. (1997) 
and Baring et al. (1999) also modeled IC 433 with synchrotron emission in the radio band and relativistic 
Bremsstrahlung in 7-rays. De Jager & Mastichiadis (1997) included inverse Compton scattering in their 
model of SNR W44. Bykov et al. (2000) have recently analyzed the non-thermal emission from a SNR 
interacting with a molecular cloud, modeling it as a highly inhomogeneous structure consisting of a forward 
shock of moderate Mach number, a cooling layer, a dense radiative shell, and an interior region filled with 
hot tenuous plasma. Particularly for SNRs with mixed morphology (remnants which are shell-like in radio 
and dominated by central emission in X-rays, Rho & Petre 1998), they found that Bremsstrahlung, syn- 
chrotron, and inverse Compton radiation of the relativistic electron population produce multiwavelength 
photon spectra in quantitative agreement with radio and high-energy observations. These are only some 
of the works devoted to high-energy emission from SNRs published in recent years. Differentiating the 
7-ray emission produced by ions from that originating in leptons is crucial for determining the origin of 
cosmic-ray nuclei (for some recent reviews and more references the reader is referred to Volk 2001, 2002; 
Drury et al. 2001; Kirk & Dendy 2001). 

After introducing a simple theoretical model for evaluating the possible hadronic 7-ray emission from 
SNRs and nearby clouds, we characterize the sample to be investigated, discuss the 7-ray flux variability of 
the sources, and study the possibility that pulsars might be possible counterparts. For each SNR-EGRET 
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Figure 4: Expected distribution of source locations (Galactic coordinates) for one year survey of the LAT 
experiment onboard GLAST. Courtesy of the GLAST Science Team and NASA. 



source pair we review and analyze the different scenarios proposed as an explanation of the 7-ray emission. 
We present CO(J=1-0) mm wavelength observations to evaluate whether there are molecular clouds in the 
vicinity of the SNRs and estimate the 7-ray flux that would be produced in each case via Tr'^-decays. Some 
new radio continuum maps are also presented. These latter maps have been processed to eliminate, as far 
as possible, the galactic contaminating diffuse emission. 

Our aim with this review is to provide a quantitative basis to analyze the possible 7-ray production 
in SNRs, providing the reader with useful information to guide future studies. Specifically, the role of 
INTEGRAL, AGILE, and GLAST satellites, and the Cerenkov telescopes HESS, VERITAS, MAGIC, and 
CANGAROO III is discussed. Several target candidates for observations with all these telescopes and satel- 
lites are mentioned. As an example, in the GeV band, the Gamma-ray Large Area Telescope (GLAST), 
which will be launched in a few years, is expected to detect ~ 10"^ high-energy 7-ray sources, thousands of 
them belonging to our Galaxy (see Figure 0]). A technical Appendix quotes the main features of GLAST, 
as well as its predecessors, AGILE and INTEGRAL, for quick reference. It is also important to clarify what 
this review is not about. Many authors have studied the evolution of SNRs or their emission properties 
from very sophisticated numerical modeling points of view, during the last years. We shall not particularly 
focus on those, except briefly when dealing with the case by case analysis for specific SNR-EGRET source 
pairs. 

The rest of this work is organized as follows. In the next Section we introduce a simple model to 
account for the hadronic 7-ray emission in SNRs and their neighborhoods. Section 3 refers to Relativistic 
Bremsstrahlung as a competing process. Section 4 analyzes the spectral changes that difussion could 
produce on the observed 7-ray spectrum. The general characteristics of the SNR sample that we shall 
analize and the possible pulsar counterparts are presented in Sections 5, 6, and 7. The variability in the 
7-ray emission for the Third EGRET sources under study is assesed in Section 8. Section 9 gives a brief 
account of the observations and data extraction techniques used in this review. A case by case analysis 
of all coincident pairs between SNRs and 7-ray sources is given in Section 10. Some particular cases in 
which SNRs were discovered by their high energy emission are discussed in Section 11. The TeV-emission 
properties and the prospects for new observations using new TeV-telescopes are discussed in Section 12. 
Finally, Section 13 presents a very brief overview and some concluding remarks. 

2 Phenomenological model for the hadronic 7-ray emission in SNRs 
and their environs 

We first present a simple model for the hadronic 7-ray emission from "bare" SNRs, and those interacting 
with molecular clouds. Further details can be obtained from Dorfi (1991, 2000), Drury et al. (1994), Aha- 
ronian et al. (1994), and Aharonian &; Atoyan (1996). We shall partially follow Morfill et al. (1984) and 
Combi & Romero (1995) in our presentation. Our intention is not to arrive at the most precise theoretical 
model for an individual SNR, but rather to have a simple, straightforward and robust, albeit crude, method 
of obtaining and inter-comparing 7-ray fluxes due to nucleus- nucleus interactions in interacting SNRs. 
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Let us consider the expansion of a SNR in a homogeneous medium. If this expansion is adiabatic, we 
can use the Sedov's solutions (1959), which give the time since the explosion and the velocity of the shock 
front, respectively, as 

t ~ 1.5 X 10\]/^ E^I^^rI^\i, (1) 

Vs ~ 21.6 X W'^nZY^El(^R^^^\m s'^ (2) 

Here, E^i is the energy of the SN explosion in units of 10^^ erg, Ri is the SNR radius in units of 10 pc, 
and n_i is the medium density in units of 0.1 cm~^. The CR energy per time unit incoming to the SNR is 

E{t) = k, ecR 47ri?s(t)' v,{t), (3) 

where the dot means derivative with respect to time, Rs is the SNR time-dependent radius, ecR is the 
background CR ambient density (~ 1 eV cm~^ in the solar neighborhood), and is the enhancement 
factor due to re-acceleration by Fermi mechanism at the shock front (see Jones 2001 for a recent review and 
references on acceleration details). If we assume equipartition, i.e. that the energy flux from the unshocked 
medium is converted in equal parts into electromagnetic energy, thermal energy, and CR enhancement 
(Morfill et al. 1984), we can write the previous expression as 

3 1 — z/inv^ 

where is the downstream to upstream ratio of kinetic energy flux in the shock frame (^ ~ 0.06 for strong 
shocks), n is the mean molecular weight, and n is the unshocked particle density. This expression states 
that the power available for accelerating CRs is 1/3 of the mechanical energy flux across the shock. 

When the SNR expands from a radius comprising a volume V{ti) to a volume V{t2), the energy 
decreases accordingly as 

Eit2) _ fv{h)y-' 

Eih) [V{t2)j ' 

with the adiabatic index being 7 = 4/3, Using the previous equations, the CR energy in the SNR between 
times ti and t2 is 

Through the Sedov solutions, this leads to a ratio 



-E'sN 5 y \t2j 

We shall adopt t2 as the actual age of the SNR, estimated from observational data and the Sedov solutions, 
and we shall assume the initial time ti as that obtained when the SNR has swept about SAfg of interstellar 
material, starting then the Sedov phase (Lozinskaya 1992, pp. 205ff). The radius at which this happens is 
i?s = (3 • 5M0/47r/Ltm/fno)^/^, and typical values for ti are in the range 200-2000 years. 

In simplified models of SNRs, the remnant is divided into three regions: an interior region filled with 
hot gas and accelerated particles but very little mass, an immediate post shock region where most of the 
matter is concentrated, and a shock precursor region where the accelerated particles diffusing ahead of the 
shock affect the ambient medium. Following Drury et al. (1994) the production rate of 7-rays per unit 
volume can be written as 

= £^n = q^nEcR, (8) 
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where n is the number density of the gas, Eqr is the CR energy density, and q-y is the 7-ray emissivity nor- 
mahzed to the CR energy density, = £^/Ecr- The total gamma-ray luminosity is given by / q^nEcRd^r, 
which can be written as q^ {MiEq^i + M2ii'cR2)) where Mi is the total mass in the precursor region, M2 
that in the immediate post-shock region, and E'en 1,2 are the corresponding CR energy densities. Since 
particle diffusion occurs across the shock front, we have i?cRi = -£'CR2- This value is also probably not 
very different from Eqr^^, the energy density in the interior of the remnant, because of two reasons. First, 
there is diffusive coupling between the acceleration region around the shock and the interior of the remnant 
(Drury et al. 1994). Second, if the acceleration is efficient, CRs provide a substantial, if not the dominant, 
part of the interior pressure and the interior of the remnant has, for dynamical reasons, to be in pressure 
equilibrium. It follows that, to order of magnitude, the CR energy density throughout the remnant and in 
the shock precursor can be taken as Eqri = -E^CR2 ~ Eqrs ~ 30i?sN /47ri?^ , where 9 is, again, the fraction 
of the total supernova explosion energy, Eq,n, converted to CR energy and R is the remnant radius. Thus, 
the 7-ray luminosity results 

= ,,(M, + M.)^ . eq.,Es,n . 10^^ (^) (^) ph s~\ (9) 

where n is the ambient density. The exact value of 6 depends on the details of the model, for which Eq. 
((T)) gives an example. For different plausible injection models, 9 is roughly constant throughout the Sedov 
phase with only a moderate dependence on external parameters such as the ambient density (Markiewicz 
et al. 1990). If the SNR is located at a distance d, the hadronic 7-ray flux is 

F(> 100MeV)s™ ~ 4.4 x 10-« (^) {^f (^) ph cm-'^s-, (10) 

where d is the distance to the remnant. Only for very high densities can the usually observed 7-ray sources 
can be due to the remnant itself. In general, the flux provided by the previous equation is far too low to 
produce a detectable EGRET source (Drury et al. 1994), at typical galactic distances. 

However, stronger emission can be produced if there are molecular clouds in the vicinity of the SNR 
where the locally accelerated protons interact with target ions, producing pions and hence enhancing the 
7-ray flux (eg. Montmerle 1979; Dorfi 1991, 2000; Aharonian et al. 1994). The expected total flux is 

Fi = ^[. n{r)q^{f)d\. (11) 



Vb 



Neglecting all possible gradients within the cloud, this equation reduces to 

Mci q^ 



= (12) 



where M^i is the mass of the cloud. In particular, we may write 

d 



F{> lOOMeV)^^"'''^ ~ 10~^Af3 ( f- ] qJ> lOOMeV) ph cm'''s~', (13) 

Vkpc/ 

where M3 is the mass of the target cloud in units of 10^ Mq, and q^ is the 7-emissivity in units of 10~^^ 
s~^ (H - atom)^^. The factor q^ will be enhanced in comparison with its normal value because of the local 
CR source. In a passive giant molecular cloud exposed to the same proton flux measured at the Earth, the 
7-ray emissivity above 100 MeV is equal to 1.53 r/ g_25(> 100 MeV)(II — atom)~^s~^, where the parameter 
r] ~ 1.5 takes into account the contribution of nuclei both in CRs and in the interstellar medium (Dermer 
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et al. 1986; Aharonian 2001). In clouds near CR accelerators, it may be much higher than this value. If 
the shape of the CR spectrum in the cloud does not differ much from that existing near the Earth, we can 
approximate 

~ ~ ks. (14) 

^7,0 ^CR,0 

Following Morfill & Tenorio Tagle (1983), we can use Eq. ^ to obtain 

The Sedov solutions for each of the SNRs considered below should be used in this latter expression to 
obtain k^, and thus F{E > 100MeV)'=^°"'^. One immediate test of energetic consistency is to check that 
-^CR = ^secR/(4/3)7ri?'^ < E'sN, for the obtained value of k^ and the assumed value of -EsN- If the previous 
inequality is not valid, one or more of the simplifying assumptions of the model are not correct for the 
particular case under analysis. 

The expected 7-ray flux in the TeV region by a SNR is (Aharonian et al. 1994) 

F,(> E) = fr 10-1° A cm'h-\ (16) 

where the factor j4 is ^ 

(lO^^erg) (kpc) (cm^)'^ ' ^ ^ 

n is the medium density, and /r is a function of the index in the differential power-law proton spectrum (T), 
equal to 0.9, 0.43, and 0.19 for r=2.1, 2.2, and 2.3, respectively. This estimate, however, usually exceeds 
that obtained when the GeV spectral index is extrapolated up to TeV energies. When that is the case, the 
extrapolated flux (with the same spectral index) will be considered a safer estimate. To extrapolate the 
GeV flux up to TeV energies we assume (Thompson et al. 1996) 

where is a constant, E is given in MeV, and Eq is a reference energy. This constant can be obtained 
simply by integrating the flux, 

/•lOGeV / K \ -r 

K^F^^I / — dE, (19) 

JlOOMeV V-C/O/ 

where Fph is the observed total flux (that quoted in the 3EG Catalog, for instance). Once K is known, 
the flux in any given energy interval Ei - E2 is just, 

F{Ei,E2)=K j dE. (20) 

We are extrapolating, then, the measured spectral index at MeV-GeV energies assuming that there is no 
spectral change at higher energies. Actually, this assumption is a simplification not compatible with TeV 
observations of several sources (see below). In any case, this extrapolation will always provide an upper 
bound to the high-energy photon flux. 
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3 Relativistic Bremsstrahlung 



Since relativistic electron Bremsstrahlung and nucleus-nucleus induced pion-decay are competing processes 
in the generation of 7-rays from molecular clouds exposed to a nearby CR accelerator, it is necessary to 
assess the relative weight of each contribution if we are to quantitatively address the question of the pos- 
sible SNR origin of nucleonic CRs. 

The 7-ray emissivity at a given energy E from relativistic Bremsstrahlung is (e.g. Longair 1994, p.267- 
269) 

in-2i 

qB{E) = KE-P m'^ s"^ GeV'S (21) 

p — 1 

where it is assumed an electron power-law distribution, Ne{E) = KE^^^. We are interested in the 7-ray 
radiation above 100 MeV, so we integrate the previous expression to obtain 

iQ— 21 i-oo 10 ■^^ 
qB{E > O.lGeV) = -n^-s K / E-^dE m"^ s"^ = :rr^n^-z K O.l^^^"^^ m"^ s~\ (22) 



P - 1 ™ Jo.lGcV {P - 1)" 

This same population of relativistic electrons will also radiate at radio wavelengths, via the synchrotron 
mechanism. The synchrotron spectrum of a power-law electron energy distribution is (Longair 1994, p. 
261) 

J{u) = 23.44 a{p) B^^^+^^Z^ K { ^-^^S x 10^^ ^ " ^ ^^^^ 



where B is the magnetic field measured in Tesla, and 



2 (p + i)r(| + |) ' ^''^ 

is a numerical coefficient depending on the spectral index. Then, the ratio between the 7-ray flux emitted 
by relativistic Bremsstrahlung and the synchrotron emission results 



_ qB{E > lOOMeV) _ F{E > lOOMeV) _ 4.3 x IQ-^^ „-(i+p)/2 (p-i)/2 1 2 -1 

R- J^^ -^d) ^-'^'^/^G ^H. Jy cm s , (25) 



where we have defined 

b{p) = 10-^(^+^^(3.2 X ioi5)(p-i)/2 _ 1)2 (26) 
and converted units to the cgs system. 

If F(E > 100 MeV) is known, estimating the right hand side of Eq. ()25() for the measured spectral 
photon index and the derived density and magnetic field, the expected value of F^ [Jy] can be obtained. 
This is the radio emission that should be observed if the 7-rays are from relativistic Bremsstrahlung. If 
the 7-ray source is not superposed with the bulk of the synchrotron radio/X-ray emission from the SNR, 
this tends to favor a nucleus-nucleus origin of the high-energy flux, rather than a electron Bremsstrahlung 
scenario. In general, at the high densities found in molecular clouds, inverse Compton scattering can be 
ruled out as the main mechanism contributing to the 7-ray emission (e.g. Gaisser et al. 1998). Above 100 
MeV, the relevant cross sections and estimates of the electron-nucleon density ratio show that relativistic 
Bremsstrahlung dominates over inverse Compton processes (see, for instance, Stecker's 1977, Figure 1 and 
2). In particular, de Jager & Mastichiadis (1997) have shown that for molecular densities above 10 cm^'^ 
7-ray fluxes above 70 MeV are dominated by Bremsstrahlung when electrons are considered (see their Fig- 
ure 4). In what follows, since we shall mainly consider high-density scenarios, relativistic Bremsstrahlung 
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will be the main alternative to nucleonic interactions in evaluating the origin of the observed 7-ray flux in 
the MeV-TeV energy range. 



4 Diffusion of CRs and 7-ray spectral evolution 

The spectrum of 7-rays generated through vr'^-decay at a source of proton density Up is 

F,iE,) = 2 r dE^, (27) 

j£;mm y/E^ -m^ 

where 

2 

m 



ET-{E,) = E, + -^, (28) 



and 



F^E^) = ATTUp / " Jp{Ep r^^%' dEp. (29) 



V 



Here, (iiT,r(-£'7r> Ep)/dET^ is the differential cross-section for the production of vr^-mesons of energy E^^ by a 
proton of energy Ep in a p — p collision. If the proton spectrum Jp{Ep) at the 7-ray production site is 

Jp{Ep) = KEf, (30) 

we can also expect a power-law spectrum at 7-rays: 

F^{E.i) E-^ . (31) 

However, the spectrum given by Eq. (|3()|) is not necessarily the same proton spectrum at the acceleration 
site. If there is diffusion, we shall have, instead 

Jp{Ep, r, t) = ^f, (32) 

where f{Ep, r, t) is the distribution function of protons at an instant t and distance r from the source. 
The distribution function satisfles the well-known diffusion equation (Ginzburg & Syrovatskii 1964): 

dl_D{Epld_,dl J_ 
dt ~ r2 d/ 5r 

where P = —dEp/dt is the continuous energy loss rate of the particles, Q = Q{Ep, r, t) is the source 
function, and D{Ep) is the diffusion coefficient, for which we assume here no dependency on r or t, i.e. the 
particles diffuse through an homogeneous, quasi-stationary medium. 

We assume that D{Ep) oc Ep and / oc E"^ with continuous injection given by Q{Ep, t) = QoE~^q{t), 
which is appropriate for a supernova remnant (Aharonian &: Atoyan 1996). Further simplicity can be 
achieved assuming that the source is constant after turning on at some instant, i.e. q{t) = for t < and 
q{t) = 1 for t > 0. Atoyan et al. (1995) have found a general solution for Eq. (|33() with arbitrary injection 
spectrum, which with the listed assumptions leads to: 
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Figure 5: Temporal and spectral evolution of CR fluxes at different (10 pc, 30 pc, and 100 pc) distances 
from an impulsive proton accelerator. A power-law proton spectrum with T = 2.2 and total energy 
Wp = 10^'^ erg are assumed. Curves 1, 2, 3, and 4 correspond to an age of the source of t = 10^ yr, 
10^ yr, 10^ yr, and 10^ yr , respectively. An energy-dependent diffusion coefficient D{E) with power-law 
index 6 = 0.5 is adopted. The left panel presents results for D = D28cm^s~^, whereas the right panel 
corresponds to D = 10^^D28 cm^s~^. The hatched curve shows the local (directly measured) flux of CR 
protons. More details are given in Aharonian & Atoyan (1996). 



In this expression, i?difr = Rdis{Ep, t) is the diffusion radius which corresponds to the radius of the 
sphere up to which the particles of energy Ep propagate during the time t after the injection. Now, for 
D(Ep) = aD2sEp, where D28 = cm^ s~^, and for -Rdiflf ^ ^) i-e- when the target is well immerse 

into the cosmic ray flux, Eq. (|3H) reduces to 



QoEp 



ir+6) 



and then, from Eq. (|32|) . we get: 



f{Ep, r) = ^^^^ , (35) 

Hence, as has been emphasized by Aharonian & Atoyan (1996), the observed 7-ray flux F^{E^) oc E-y ^^'^^'^ 
can have a significantly different spectrum from that expected from the particle population at the source 
(the SNR). Standard diffusion coefficients 5 ~ 0.3 — 0.6 can explain 7-ray spectra as steep as F ~ 2.3 — 2.6 
in sources with particles accelerated to a power-law Jp{Ep) oc E~'^ if the target is illuminated by the 7r°- 
decays are at sufficient distance from the accelerator. This can explain observed discrepancies in the particle 
spectral indices inferred from SNR at different frequencies, even if all particles, leptons and hadrons, are 
accelerated to the same power-law in the source. 



CRs with total energy Wp and injected in the interstellar medium by some local accelerator reach a 
radius R{t) at instant t. Their mean energy density is Wp ^ 0.5(VFp/10^°erg)(R/100pc)-3 eVcm'^. Thus, 
in regions up to 100 pc around a CR accelerator with Wp ~ lO^'^ erg, the density of relativistic particles 
may significantly exceed the average level of the "sea" of galactic CRs, wgcr ~ leV/cm^. In Figure El the 
differential flux of protons at distances R =10, 30, and 100 pc from an "impulsive" accelerator, with total 
energy Wp = W^^ erg are shown. The spectrum of CRs at the given time and spatial location can differ 
significantly from the source spectrum. The diffusion coefficient in this figure is assumed in a power-law 
form, D{E) oc E^'^ above 10 GeV, and constant below 10 GeV. The commonly adopted value at 10 GeV 
is about lO^^cm^s"^, however smaller values, e.g. lO^^cm^s"^, cannot be excluded, especially in active 
star forming regions (Aharonian 2001). The existence of massive gas targets like molecular clouds in these 
regions may result in 7-ray fiuxes detectable by EGRET, if (W50 • M^)/dl^^ >~ 0.1, where M5 is the 
mass of the cloud in units of IO^Mq (Aharonian 2001). In the case of energy-dependent propagation of 
CRs, large variety of 7-ray spectra is then expected, depending on the age of the accelerator, duration of 
injection, the diffusion coefficient, and the location of the cloud with respect to the accelerator. 

The comparison of 7-ray fluxes from clouds located at different distances from an accelerator may 
provide unique information about the CR diffusion coefficient D{E). Similar information may be obtained 
from a single 7-ray emitting cloud, but in different energy domains. For the energy-dependent propagation 
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Figure 6: Gamma-ray emissivities in terms of i?^ x qy(E^) at different times t and distances R from a 
proton accelerator. The right hand side axes shows the 7-ray fluxes, x F^{E^), which are expected from 
a cloud with parameter M^/d'^^^ = 1. The thin and bold curves correspond to times t = 10^ and 10^ yr, 
respectively. Fluxes at distances R = W and 30 pc are shown by solid and dashed lines. The power-law 
index, the total energy of protons, and the diffusion coefficient are the same as in the right panel of Figure 
121 The curve shown by full dots corresponds to the 7-ray emissivity (and flux) calculated for local CR 
protons. In order to take into account the contribution of nuclei, all curves should be increased by a factor 
of r/ ~ 1.5. More details can be found in Aharonian & Atoyan (1996). 



of CRs the probability for simultaneous detection of a cloud in GeV and TeV 7-rays is not very high, because 
the maximum fluxes at these energies are reached at different epochs (see Figure EI). The higher energy 
particles propagate faster and reach the cloud earlier, therefore the maximum of GeV 7-radiation appears 
at the epoch when the maximum of the TeV 7-ray flux is already over. In the case of energy-independent 
propagation (e.g. due to strong convection) the ratio F^f{> 100MeV)/F^(> 100 GeV) is independent of 
time, therefore the clouds that are visible for EGRET at GeV energies would be detectable also at higher 
energies, provided that the spectral index of the accelerated protons be F < 2.3. 

Summing up, special care must be taken in analyzing the a priori expectations for the SNR-cloudy 
medium scenario: when particles diffuse into the ISM before reaching the cloud, a quite different spectrum 
from that of the injecting source can finally emerge, depending on the distance to the cloud and the 
diffusion index. Even when molecular clouds are overtaken by the expanding shell of the SNR and there 
is clear evidence of interaction, a strong magnetic field could produce differences in the spectra. 

5 Sample and correlation analysis 

Tabled shows those 3EG sources that are positionally coincident with SNRs listed in the latest version of 
Green's Catalog (2000). From left to right, columns are for the 7-ray source name, the measured flux in the 
summed EGRET phases P1234 (in units of 10~^ ph cm~^ s~^), the photon spectral index F, the EGRET 
class of source (em for possibly extended and C for confused), the variability indices / (as in Torres et 
al. 2001a,c) and t^^^^^^^^I^ (as in Tompkins 1999), information about coincidences with radio pulsars ("y" 
stands for a pulsar within the error box), the SNR identification (including other usual names when avail- 
able), the angular distance between the center of the 7-ray source position and the center of the remnant 
(in degrees), the size of the remnant (in arcmin), and finally the SNR type T (S for shell- like emission, F 
for filled-centre or plerionic remnant, and C for composite). A separate section below analyzes the possible 
pulsar associations. All remnants were considered as circles with a radius equal to the major axis of the 
ellipse that better fits their shape, when such is given in Green (2000). It is interesting to see that many of 
the 3EG sources involved in the associations are classified as extended, and all of them as confused. Also, 
it was already noted (Romero et al. 1999a) that in the 3EG catalog not all the positional coincidences 
with SNRs are SNOBs (SNRs in OB associations), as it was the case in the studies by Montmerle (1979) 
and Yadigaroglu & Romani (1997) using previous samples. 

The adoption of the 2000 edition^ of Green's Catalog (2000) does not produce any substantial statisti- 

new edition -2001- of the catalog has recently appeared. The results presented herein are not changed when the 2001 
edition is considered. 
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Table 1: Positional coincidences between supernova remnants quoted in Green's Catalog (2000), and unidentified 3EG EGRET sources. 
See text for the meaning of the different columns. 



7-source 




r 






Class 


/ 




r 


p? 


SNR 


Other name 


A9 


Size 


T 


0542+2610 


14.7+3.2 


2, 


.67+0, 


.22 


em 


c 


3 


.16 


0.701:12 




G180.0-1.7 




2.04 


180 


S 


0617+2238i'2 


51.4+3.5 


2, 


.01+0. 


.06 


C 




1 


.68 


0.260Ji 




G189.1+3.0 


IC443 


0.11 


45 


s 


0631+0642i'3 


14.3+3.4 


2, 


.06+0. 


.15 


c 




1, 


.52 






G205. 5+0.5 


Monoceros 


1.97 


220 


s 


0634+0521 


15.0+3.5 


2, 


.03+0. 


.26 


em 


c 


1 


.02 


72.0^15 




G205.5+0.5 


Monoceros 


2.03 


220 


s 


1013-5915 


33.4+6.0 


2, 


.32+0. 


.13 


em 


c 


1 


.63 


0.22O:^« 


y 


G284.3-1.8 


MSH 10-53 


0.65 


24 


s 


1102-6103 


32.5+6.2 


2, 


.47+0. 


.21 


C 




1 


.86 






G290.1-0.8 


MSH 11-61A 


0.12 


19 


s 






















G289.7-0.3 




0.75 


18 


s 


1410-6147^ 


64.2+8.8 


2, 


.12+0. 


.14 


C 




1, 


.22 


n oo0.5,5 


y 


G312.4-0.4 




0.23 


38 


s 


1639-4702 


53.2+8.7 


2, 


.50+0. 


.18 


em 


c 


1 


.95 


U-UUo.oO 


y 


G337.8-0.1 
G338.1+0.4 
G338.3+0.0 


Kes 41 


0.07 
0.65 
0.57 


9 

15 
8 


s 
s 
s 


1714-3857 


43.6+6.5 


2, 


.30+0. 


.20 


em 


c 


2, 


.17 


n 1 50.38 


y 


G348.5+0.0 
G348.5+0.1 
G347.3-0.5 


CTB 37A 


0.47 
0.50 
0.85 


10 
15 
65 


s 
s 
s 


1734-32325 


40.3+6.7 








C 




2, 


.90 


n nno-24 

".""0.00 




G355. 6+0.0 




0.16 


8 


s 


1744-3011 


63.9+7.1 


2, 


.17+0. 


.08 


C 




1 


.80 


0.38[j:^o 




G359.0-0.9 




0.41 


23 


s 






















LxoOy.i — U.O 




u.zo 


OA 


Q 




1746-2851** 


119.9+7.4 


1, 


.70+0. 


.07 


em 


c 


2 


.00 


o.5o8:i 




GO.0+0.0 

GO. 3+0.0 




0.12 

0.19 


3.5 

16 


s 
s 


1800-2338i'^ 


61.3+6.7 


2, 


.10+0. 


.10 


C 




1 


.60 




y 


G6.4-0.1 


W28 


0.17 


42 


c 


1824-1514 


35.2+6.5 


2, 


,19+0. 


.18 


C 




3 


.00 


".""0.00 


y 


G16. 8-1.1 




0.43 


30 




1837-0423 


<19.1 


2, 


.71+0. 


,44 


C 




5, 


.41 


12.0fi7 


y 


G27.8+0.6 




0.58 


50 


F 


1856+0114« 


67.5+8.6 


1, 


,93+0, 


,10 


em 


c 


2, 


.92 


o.soJiiJ 


y 


G34.7-0.4 


W44 


0.17 


35 


s 


1903+0550'' 


62.1+8.9 


2, 


,38+0, 


,17 


cm 


c 


2, 


.28 


0.35°:?° 


y 


G39.2-0.3 


3C396, HC24 


0.41 


8 


s 


2016+3657 


34.7+5.7 


2, 


,09+0, 


,11 


C 




2, 


.06 


A 07O.75 




G74.9+1.2 


CTB 87 


0.26 


8 


F 


2020+4017^'^ 


123.7+6.7 


2, 


.08+0. 


.04 


c 




1, 


.12 


U-U<o.OO 


? 


G78.2+2.1 


W66, 7-Cygni 


0.15 


60 


s 



1 Association proposed by Sturner & Dermer (1995) and Esposito et al. (1996). ^ GeV J0617+2237 ^ GeV J0633+0645. ^ Association proposed by 
Stumer & Dermer (1995). ^ GeV J1732-3130. ^ GeV J1746-2854. ^ GcV J1800-2328. ^ GcV ,11856-0115. ^ GeV J2020+4023. GeV sources 

compiled in the GeV ASCA Catalog (Roberts et al. 2001a). 



Table 2: Evolution of the number of positional coincidences between SNRs and unidentified EGRET 
sources. In the last row, 6 hkely artifacts are disregarded in the 3EG Catalog. 



Catalog Unidentified Real Number of SNRs Significance 

EGRET detections coincidences in Green's catalog (statistical) 



First EGRET Catalog^ 


37 


13 (35%) 


182 


1.8cr 


2EG'' 


32 


7 (22%) 


194 


? 


2EG" 


32 


5 (16%) 


14S 


? 


2EG'^ 


33 


10 (30%) 


194 


h 


SEG'' 


81 


22 (27%) 


220 


5.7(7 


3EGf 


75 


19 (25%) 


220 


4.8(7 



a: Sturner & Dermer (1995). b: Sturner, Dermer & Mattox (1996). c: Esposito, Hunter, Kanbach & Sreekumar 

(4-996). d: Yadigarpglu & Romani X1997),, e: Romero, Benaglia & Torres (1999)„f: Torres^t,al„(2G01b). g:. Only 
l^ieure 7: Left."^ Statistical results lor the random association between SNKs and EGKEl sources at low- 
raaio-hright^NiRs, ilux at 1 ,GHz greater than lOO Jy, werCpConsidcrcd, h: Comnuted for pairs. , „ 
latitudes. Kighf: Distribution of the 7-ray fluxes as a function of tne photon specfTral index. Iwo sources 

seem to differentiate from the rest. One of these sources (see below) has been recently identified with a 

blazar. 

cal difference with respect to the previous edition of 1998. Only 5 SNRs were added. However, there is a 
notable particular difference in the case of the EGRET source 3EG J1714-3857, which now coincides with 
three supernova remnants instead of two. One of these SNRs (the new one in Green's catalog, G347.3-0.5) 
appears to be amongst the strongest cases for SNR shock/7-ray/nucleonic cosmic-ray source associations 
known to date (Butt et al. 2001, see below). In addition, 2EG J1801-2312 shifted its position half a degree 
when converting into 3EG J1800-2338, also affecting previous positional coincidences. The evolution in 
the number of coincidences between SNRs and EGRET sources since the First EGRET Catalog until the 
current situation is shown in Table El (Torres et al. 2001b). 

The Poisson probability for the 19 coincidences to be a chance effect is 1.05 x 10~^, i.e., there is an 
a priori 0.99998 probability that at least one of the positional associations in Table ^ is physical. This 
expected chance association was computed using thousands of simulated sets of EGRET sources, by means 
of a numerical code described elsewhere (Romero et al. 1999a, b, Sigl et al. 2001). Figure [7| shows the result 
of numerical simulations for these random populations of 7-ray sources. Figure |H1 shows the distribution 
of the 7-ray photon spectral index for the sample of 3EG sources coincident with SNRs. Some cases of 
possible physical associations mentioned in the literature and discussed in the next sections are indicated. 

If some SNRs interact, as expected, with nearby massive clouds producing enhanced 7-ray emission 
through hadronic/Bremsstrahlung interactions, cases in which there is just a marginal coincidence between 
the center-points of the SNRs and the centers of the EGRET sources should be also considered, since the 
peak 7-ray emissivity will likely be biased towards the adjacent cloud. So we have also looked at the 
positional coincidences between unidentified EGRET sources and the region just around the SNRs. We 
have done so by artificially enlarging the size of the SNR by half a degree. We found that there are 26 
coincidences of this kind, including those in Tabled Then, there are 7 new cases perhaps worthy of further 
study. These cases are shown in Table El Interestingly, the expected chance coincidence in this case is at 
the level of 14.8ib3.14, still Aa lower than the real result, implying a probability of 2.5 x 10^^ for the real 
result to be a random (Poisson) fluctuation. In order to quantify the role played by the 7 new sources 
in this result new simulations were carried out, now considering only these sources. The chance result is 
2.3ibl.2, again 4fT below the real number of coincidences. In the present review, however, we shall focus 
only on those SNRs that present positional correlation with 7-ray sources, i.e. only the cases listed in 
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Figure 8: Distribution of the 7-ray photon spectral index. Some SNR-3EG coincidences for which the 
physical association has been suggested in the literature are indicated. 



Table 3: Supernova remnants and nearby (but not coincident) unidentified 3EG. 



7-source T Class / t SNR Size T Other 



0229+6151 


39.9+6.2 


2.29+0.18 


C 




1.3 


•?70-74 


G132.7+1.3 


1.50 


82 


S 


Of/OB 


1736-2908 


51.5+9.1 


2.18+0.12 


C 




2.4 


U.DDo.40 


G359. 1+0.9 


0.73 


12 


S 




1741-2050 


24.1+3.9 


2.25+0.12 


C 




2.1 


410-70 


G6.4+4.0 


1.00 


31 


s 




1823-1314'''^ 


102.6+12.5 


2.69+0.19 


c 




2.9 


n 79I.37 


G18. 8+0.3 


0.88 


17 


s 


OB 


1826-1302 


66.7+10.1 


2.00+0.11 


c 




2.6 


7^1-28 


G18.8+0.3 


0.88 


17 


s 


OB 


1928+1733 


157.0+36.9 


2.23+0.32 


em 


C 


3.9 


0.82q;43 


G54.1+0.3 


1.21 


1.5 


F? 




1958+2909= 


26.9+4.8 


1.85+0.20 


em 


C 


1.6 


4'^0.98 
"•^'-'0.15 


G65. 1+0.6 


1.36 


90 


S 





a : Association proposed by Sturner & Dermer (1995), and Esposito et al. (1996). b : GeV J1825— 1310. c : GeV 
J1957-8859. 

Tabled 



6 SNRs coincident with 7-ray sources 

We now analyze the sample of SNRs in Table d in more detail. Table |1] presents the radio fluxes of the 
SNRs, together with their spectral index and known distance estimates (with the corresponding refer- 
ences). Distances are only approximate since several different values for the same SNR can be found in 
the literature. When no direct determination is available, estimates can be made using the radio surface 
brightness-to-diameter relationship, known as T, — D (Clark & Caswell 1976, Milne 1979, Case & Bhat- 
tacharya 1998). For these cases, marked with a star in Table |3 the distances given by the new S — D 
introduced by Case & Bhattacharya (1998) will -unless otherwise noted- be adopted. A double star symbol 
means that neither a distance determination nor an estimate is available. There is only one such case in 
Table 01 for which the distance to a coincident OB association was assumed. Thus, distances in Tabled 
marked with one or two stars are, respectively, more uncertain than the others. 

Using the estimated distance to each remnant in Table 0J we have calculated the approximate intrinsic 
7-ray luminosity of the putative region producing the high-energy source in the energy range 100 MeV-10 
GeV using the observed EGRET flux and photon indices (see Table d), assuming isotropic emission. We 
have also re-ordered Green's Catalog according to descending radio flux, and the rank in this list is given 
for each remnant. It is interesting that from the first 20 SNRs with the highest radio fluxes, only 6 appear 
to be correlated with EGRET sources. Sturner &: Dermer (1996) have noted that SNRs not correlated with 
2EG sources were either more distant than 7-Cygni (G78. 2+2.1) and IC 433 (G189. 1+3.0), or presenting 
a far smaller radio flux. This trend is not observed now with the larger 3EG sample. 
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Table 4: Properties of the SNRs coincident with 3EG sources. The distances quoted are reported or 
discussed in the cited references. Radio fluxes and spectral indices (a, such that oc i/") are taken from 
Green's (2000) Catalog. 



7-source 




a 

kpc 


Ker. 


radio 

Jy 


a 


HanK; 


T 

J-J-y 

erg e 


-1 


0542+2610 


G180.0-1.7 


0.8-1.6 


1 


65 


varies 


24 


9.65 


1033 


0617+2238 


G189.1+3.0 


1.5 


2 


160 


0.36 


11 


1.01 


1035 


0631+0642 


G205.5+0.5 


0.8-1.6 


3 


160 


0.5 


12 


1.70 


1034 


0634+0521 


G205.5+0.5 


0.8-1.6 


3 


160 


0.5 


12 


1.85 


1034 


1013-5915 


G284.3-1.8 


2.9 


4 


11 


0.3? 


91 


1.71 


1035 


1102-6103 


G290. 1-0.8 


7 


5 


42 


0.4 


34 


8.46 


1Q35 




G289.7-0.3 


8.2 


* 


6.2 


0.2? 


123 


1.11 


1036 


1410-6147 


G312.4-0.4 


1.9-3.1 


6 


45 


0.26 


32 


3.06 


1035 


1639-4702 


G337.8-0.1 


12.3 


7 


18 


0.5 


67 


4.17 


1036 




G338.1+0.4 


9.9 


* 


4? 


0.4 


153 


2.70 


1036 




G338.3+0.0 


8.6 


* 


7? 


? 


121 


2.04 


1036 


1714-3857 


G348. 5+0.0 


11.3 


8 


9 


0.4? 


96 


3.47 


1036 




G348.5+0.1 


11.3 


8 


72 


0.3 


21 


3.47 


1036 




G347.3-0.5 


6.3 


9 


? 


? 


? 


1.07 


1036 


1734-3232 


G355.6+0.0 


12.6 


* 


3? 


? 


170 






1744-3011 


G359.0-0.9 


6 


10 


23 


0.5 


54 


1.65 


1036 




G359.1-0.5 


8.5-9.2 


8-11 


14 


0.4? 


81 


3.56 


1036 


1746-2851 


GO.0+0.0 


8.5 


8 


100? 


0.8? 


18 


1.20 


1037 




GO.3+0.0 


8.5 


12 


22 


0.6 


55 


1.20 


1037 


1800-2338 


G6.4-0.1 


1.6-4.2 


13 


310 


varies 


6 


4.04 


1035 


1824-1514 


G16.8-1.1 


1.48 




2? 


? 


186 


5.42 


1034 


1837-0423 


G27.8+0.6 


2 


14 


30 


varies 


44 


3.40 


1034 


1856+0114 


G34.7-0.4 


2.5 


15 


230 


0.30 


7 


4.14 


1035 


1903+0550 


G39.2-0.3 


7.7-9.6 


15 


18 


0.6 


68 


2.12 


1036 


2016+3657 


G74.9+1.2 


10 


15 


9 


varies 


104 


2.75 


1036 


2020+4017 


G78.2+2.1 


1.7 


16 


340 


0.5 


5 


3.20 


1035 



1. Anderson et al. (1996) 2. Fescn (1984) 3. Jaffc ct al. (1997) and Hcnsbcrgc ct al. (2000) 4. Ruiz & May (1986) 
5. Kaspi et al. (1997) 6. Caswell & Barnes (1985), Case & Bhattacharya (1999) 7. Koralesky et al. (1998) 8. Green 
et al. (1997), see also Reynoso & Mangum (2000) 9. Slane et al. (1999) 10. Bamba et al. (2000) 11. Uchida, Morris 
& Yusef-Zadeh (1992) 12. Kassim & Frail (1996) 13. Frail et al. (1993) and Clark & Caswell (1976) 14. Reich, 
Furst & Sofnc (1984) 15. Green (2000) and Caswell ct al. (1975) 16. Lozinskaya ct al. (2000) * From the S - D 
relationship presented by Case & Bhattacharya (1998) ** Distance assumed equal to a coincident OB association, 
Romero, Benaglia & Torres (1999). 
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Table 5: Properties of the 7-ray pulsars detected by EGRET. Pulsar parameters and distances are taken 
from Kaspi et al. (2000), except PSR B1055— 52, for which a smaller value of distance was also considered 
(Ogelman & Finley 1993, Combi et al. 1997), and Vela (Caraveo et al. 2001 and references therein), 
r = P/2P, and E = A-k'^IP/P'^, with / = lO^^gcm^. The "P1234" 7-ray fluxes and spectral indices are 
from the 3EG catalog (Hartman et al. 1999). 



Pulsar/3EG source 


P 


r 


E 


d 


F3EG[x10-8] 




V 




(ms) 


(kyr) 


(ergs~i) 


(kpc) 


(ph cm^^ s^^) 




(lOOMeV- 
-lOGoV) 


Crab/0534-2200 


33 


1.2 


5.0 X 10^** 


2.0 


226.2±4.7 


2.19±0.02 


0.01% 


Vela/0834-4511 


89 


12.5 


6.3 X 10^6 


0.25 


834.3±11.2 


1.69±0.01 


0.08% 


B1951+32/- 


39 


100.0 


3.7 X 10^6 


2.4 






0.3% 


B1706-44/1710-4439 


102 


15.8 


3.1 X 10^6 


1.8 


111.2±6.2 


1.86±0.04 


1% 


Gcminga/G633+1751 


237 


316.2 


3.1 X 10^4 


0.16 


352.9±5.7 


1.66±0.01 


3% 


B1055-52/1058-5234 


197 


501.1 


3.1 X 10^4 


0.5/1.5 


33.3±3.8 


1.94±0.10 


2/19% 



7 Pulsars within the EGRET error boxes 

Since both molecular clouds and pulsars can produce 7-rays, and because both often lie close to SNRs, it 
is important to explore the possible origin of the high-energy flux in neighboring pulsars. 

In order to firmly identify a pulsar as the origin of the 7-rays from an EGRET source, 7-ray pulsations 
must be detected at the pulsar period. However, this is not always possible because of the usually low pho- 
ton counts observed in most cases. There are only 6 high-energy 7-ray pulsars that are already confirmed, 
5 of them have associated 7-ray sources in the 3EG Catalog. We give their properties in Table |S1 Other 
candidates are usually judged by comparison with the properties of the known EGRET pulsars. 

The results of a correlation analysis between the 3EG sources superposed to SNRs (Table 
^) and pulsars are presented in this Section. The latter were extracted from the Princeton 
Catalog (Taylor et al. 1993, available on line at |http://pulsar.princeton.edu/ftp/pub/catalog/ 1 
and from the recently -partially- released Parkes Multibeam Survey (Manchester et al. 2001, 
,http : / / www . at nf . C S IRO . AU / research /pulsar/ pmsurv / ) . Adding up both surveys, there are more than 
1000 known pulsars. 

Table 13 presents the results of the correlative 3EG-radio pulsar spatial coincidence analysis: name of 
the 3EG source, name of the pulsar found within the error box, their angular separation, and the size of 
the 95% confidence contour of the 3EG source. Apart from the results quoted in Table [71 recently, Torres 
and Nuza (2003) discovered five new coincidences between EGRET sources and pulsars in the 2003 version 
of the Parkes Catalog, but they discarded their possible association based on simple spin-down energetics. 
We refer the reader to their paper for further details.^ We provide also the available information about 
the pulsar: its galactic coordinates, distance, characteristic time r = P/2P, with P and P the period and 
period derivative, respectively, spin-down energy release E = Att^IP/P^ (assuming a neutron star moment 
of inertia / = lO'^^gcm^), and the efficiency in converting the spin-down luminosity into 7-rays, if the 
pulsar alone were responsible for generating the 3EG source fiux. 



^Discussion on how gamma-ray observations of Parkes' pulsars can help distinguishing the outer gap from other models for 
gamma-ray emission can be found in Torres and Nuza (2003) and references cited therein. 
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The 7-ray efficiency was estimated as 

r? = L^/E = f^TiSF-^lE, (37) 

where is the observed 7-ray flux between 100 MeV and 10 GeV, and / is the 7-ray beaming fraction 
(0 < / < 1). This fraction is essentially unknown (see, e.g., Yadigaroglu & Romani 1995, Romero 1998), 
but it is common practice to assume / = l/47r (e.g. Thompson 2001, Kaspi et al. 2000, D'Amico et al. 
2001, Torres et al. 2001d). These efficiencies are uncertain also because they suffer (quadratically) the 
imprecise knowledge of the distance to the pulsar. 

For the confirmed 7-ray pulsars, the observed efficiencies are in the range t] €z 0.01%, ~ 3-19%), 
where the interval of upper limits is caused by the different estimates of the distance to PSR B1055— 52 
(see Ogelman & Finley 1993, Combi et al. 1997, Romero 1998). It could be considered reasonable that a 
pulsar generates 7-rays with efficiencies in the range ry € (~ 0.01%, ~ 10%). For higher values the pulsar 
would be too close to the so-called death line, where the high-energy emission is quenched (Usov 1994). 

Two of the pulsars in Table 0, PSR J1800-2343 and PSR J1902+0615, lack a confident determination 
of the period Derivative; for these it is impossible to assess their expected efficiencies. Based on the re- 
quired efficiencies and derived spin-down luminosities of the other members of the group shown in the first 
(Princeton) panel, none of the positional associations seem likely except that between 3EG 1856+0114 and 
PSR J1856+0113. For this case, the pulsar should have an efficiency of 13% in converting rotational energy 
into 7-rays. PSR J1856+0113 was already mentioned as possibly associated with the EGRET source by 
de Jager & Mastichiadis (1996). 

The bottom panel of Table (jTj) contains data that were analyzed elsewhere (Torres et al. 2001d, D'Amico 
et al. 2001, Camilo et al. 2001); here are the main conclusions. 1. The physical association between 3EG 
J1013-5915 and PSR J1016-5857 is possible, as ffist noted by Camilo et al. (2001). 2. The source 3EG 
J1410-6147 and either PSR J1412-6145 or J1413-6141 might be physically associated if either of the latter 
lie closer than their dispersion measure distances, say at the estimated distance to G312.4-0.4 (2-3 kpc). 3. 
The pulsar PSR J1637-4642 could contribute part of the high-energy budget of the source 3EG J1639-4702. 
In addition to this, the pulsars J1015-5719, J1420-6048 and J1837-0604 are also possible counterparts for 
their respective EGRET sources. 

8 Variability 

Ift SNRs, molecular clouds, or pulsars are responsible for some EGRET sources, we would expect them 
to be non-variable on the time scale of EGRET observations, i.e. from weeks to a few years. Hence, 
variability analysis of the 7-ray emission is an important tool to test the original hypothesis, in the sense 
that variable sources could be ruled out as being produced by SNRs (or pulsars). 

Three variability indices for EGRET sources have been introduced in the literature so far. The ffi'st of 
them, dubbed V, was presented by MacLauglin et al. (1996), who computed it for the sources contained in 
the Second EGRET Catalog. This method was later used, also, for short timescale studies by Wallace et al. 
(2000, 2002). The basic idea behind V is to find from the measured fluxes, and to compute V = — log Q, 
where Q is the probability of obtaining such a if the source were constant. Several critiques have been 
mentioned concerning this procedure, among them, that the scheme gets complicated when the fluxes are 
just upper limit detections. It can be shown that sources which have upper limits included in the analysis 
will have a lower V than what is implied by the data (Tompkins 1999). In addition, a source can have 
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Table 6: Positional coincidences between 3EG unidentified sources superposed to SNRs and pulsars in tlie 
Princeton Catalog and in the recently -partially- released Parkes Multibeam Survey. We show measured 
and derived pulsar parameters as well. See text for details. In the case of 3EG J1410-6147, second and 
third values of efficiency are given for two different estimates of the distance to G312.4-0.4, 1.9 kpc and 3.1 
kpc, respectively. All efficiencies are based on the new data of the 3EG Catalog (Hartman et al. 1999). 



3EG J 


PSR J 


A6» 


e 


{l,h) 


d 


T 


P 


P 


E 


V (%) 






deg 


deg 


deg 


kpc 


kyr 


ms 


10-15 


erg 


(beamed) 


Princeton 


1013-5915 


1012-5857 


0.29 


0.72 


283.7,-2.1 


10.1 


741 


820 


17.69 


1.2 xlO^^ 


> 100 


1639-4702 


1640-4715 


0.29 


0.56 


337.7,-0.4 


7.2 


188 


518 


42.02 


1.3 xlO^^ 


> 100 


1800-2338 


1800-2343 


0.12 


0.32 


6.1,-0.1 


4.8 




1030 








1824-1514 


1825-1446 


0.46 


0.52 


16.8,-1.0 


5.4 


195 


279 


22.68 


4.1 xl034 


> 100 


1837-0423 


1836-0436 


0.28 


0.52 


27.1,1.1 


4.6 


33 


354 


1.66 


1.5 xlO^^ 


> 100 


1856+0114 


1856+0113 


0.05 


0.19 


34.5,-0.5 


3.3 


10 


267 


208.408 


4.3 xlO^s 


13 


1903+0550 


1902+0556 


0.26 


0.64 


39.5,0.2 


3.9 


912 


746 


12.896 


1.2 xlO^^ 


> 100 




1902+0615 


0.48 




39.8,0.3 


10.1 




673 








Parkes 


1013-5915 


1016-5857 


0.16 


0.72 


284.1,-1.9 


3.0 


21 


107 


0.806 


2.6 X 10^** 


0.5 




1013-5934 


0.31 




284.1,-2.6 


11.3 


12561 


442 


0.557 


2.5 X 10^2 


>100 


1014-5705 


1015-5719 


0.30 


0.67 


283.69,-0.58 


4.9 


39 


140 


57.4 


8.2 xlO^s 


5 


1410-6147 


1412-6145 


0.15 


0.36 


312.3,-0.3 


9.3 


50 


315 


98.7 


1.2 X 10^5 


> 100/12/30 


1420-6038 


1420-6048 


0.16 


0.33 


313.54,+0.23 


8 


13 


68.2 


82.85 


1.0x10^^ 


2 




1413-6141 


0.28 


0.36 


312.4,-0.3 


11.0 


13 


286 


333.4 


5.7 X 10^5 


80/3/6 


1639-4702 


1637-4642 


0.46 


0.56 


337.8,+0.3 


5.8 


41 


154 


59.2 


6.4 X 10^5 


12 




1640-4648 


0.37 




338.1,-0.2 


6.1 


3501 


178 


0.806 


5.6 X 10^^ 


»100 




1637-4721 


0.45 




337.3,-0.1 


5.9 


4160 


1165 


4.44 


1.1 X 10^2 


> 100 


1714-3857 


1713-3844 


0.30 


0.51 


348.1,0.2 


6.5 


143 


1600 


177.41 


1.7 X 10^3 


> 100 




1715-3903 


0.23 




348.1,-0.3 


4.8 


117 


278 


37.688 


6.9 X 10^^ 


72 


1837-0423 


1838-0453 


0.50 


0.52 


27.1, +0.7 


8.2 


52 


381 


115.7 


8.3 X 10^4 


< 55 


1837-0606 


1837-0559 


0.14 


0.19 


26.0,+0.38 


5.0 


1045 


201.0 


3.304 


1.5x10^4 


> 100 




1837-0604 


0.16 




25.96,+0.27 


6.2 


34 


96.3 


45.170 


2.0 xlO^'^ 
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a large value of V because of intrinsic reasons or because of small error bars in the flux measurements. 
Similarly, a small value of V can imply a constant flux or big error bars. Each value of V is obtained 
disregarding those of a control population. Then, there could be pulsars with very high values of F, or 
variable AGNs with very low ones. Hence, the tisc of V to classify the variability of 7-ray sources seems not 
to produce confident results. Two other indices have been computed for all 7-ray sources: the / index and 
the r index (Torres et al. 2001a; Tompkins 1999, respectively). See Torres et al. (2001c) for a comparison 
among the results obtained with them. The idea behind the index / is to carry out a direct comparison of 
the flux variation of any given source with that shown by pulsars, which is considered as instrumental. It 
basically establishes how variable a source is with respect to the pulsar population. Contrary to Tompkins' 
index, the /-scheme uses only the publicly available data of the 3EG Catalog, and is defined as follows. 
Firstly, a mean weighted value for the EGRET flux is computed as 
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Here A^vp is the number of single viewing periods for each 7-ray source, F{i) is the observed flux in the 
i*'^-pcriod, whereas e{i) is the corresponding error in the observed flux. These data arc taken directly from 
the 3EG Catalog. A fluctuation index, n, is defined as (e.g. Romero et al. 1994): /j, = 100 x o-gd x {F)~^ , 
where asd is the standard deviation of the flux measurements. This fluctuation index is also computed for 
the confirmed 7-ray pulsars in the 3EG Catalog, and then the averaged statistical index of variability, /, 
is introduced by / = /^source 

/ < jU >puisars- We refer the reader to the work by Torres et al. (2001a,c) for 
details on the error estimates. 

Since the /-scheme is a relative classification, a result like / = 3.0 says that the flux evolution is 
three times more variable than (equivalently, 4a above) the mean flux evolution for pulsars. Torres et al. 
(2001c) mentioned that in order to get more reliable results under the /-scheme at low-latitudes it seems 
safer to consider a restrictive criterion. For instance, a source will be considered very variable only when 
I — SI> [Ip] + 3apj. Here 61 ~ 0.5/, [Ip] is the mean value of / for pulsars, apj is their standard deviation, 
and [Ip] +3(Jpj = 2.5. Rephrasing the previous constraints just in terms of /, a source will be very variable 
if / > 5. This represents a deviation of 8a from the mean /-value for pulsars. With this criterion, 3EG 
J1837-0423, with / = 5.41, is a very variable source. We have 7 non- variable sources, those having / < 1.7 
in Table 1. The rest of the sources listed there are dubious under this classification scheme. 

Tompkins (1999) used the 145 marginal sources that were detected but not included in the final official 
3EG list, and, simultaneously, all the detections within 25 deg of the source of interest. The maximum 
likelihood set of source fluxes was then re-computed. From these fluxes, a new statistic measuring the 
variability was defined as r = a/n, where a is the standard deviation of the fluxes and fi their average 
value. The final result of Tompkins' analysis is a table listing the name of the EGRET source and three 
values for r: a mean, a lower, and an upper limit (68% error bars). 

The mean value of r for pulsars -again an assumed non-variable population- is very low, 0.09, but the 
mean of the upper limits is ~ 0.2. So pulsars are consistent with having values of r up to 0.2. The deviation 
for the mean value of pulsars is 0.08. A source will be likely variable under the r scheme when the lower 
limit is at least 0.6, 3a above the mean value of the r upper limit for pulsars. A source will be considered 
non-variable when the upper limit for r is below that threshold. Sources not fulfilling either classification 
will be considered as dubious. This is also consistent with the fact that the mean r value for the known 
population of AGNs is 0.9. Using this scheme, too, 3EG J1837-0423 is a variable source. In addition, 
3EG J0631-F0642 and 3EG J0634-F0521 are also variable under r. However, Tompkins (1999) adds a word 
of caution for these two sources: the fitted flux for them is zero. Such sources arc most likely variable, 
but unknown instrument systematics or numerical problems within the r scheme could conceivable change 
these results. Also important is to note that many sources have a dubious classification: within the 68% 
error bars on r, they can be as variable as an AGN, or as non-variable as a pulsar. This is, unfortunately, 
a common situation for many of the sources. For those ones, in particular, the /-index scheme can provide 
some additional information. 

Based on the variability of the 7-ray flux, then, it is highly unlikely that 3EG J 1837-0423 is caused by 
SNR G27.8-f-0.6, or by the pulsars PSR J1836-0436 (from Princeton Catalog) and PSR J1838-0453 (from 
Parkes Catalog). In addition, the spectral index for this 3EG source is very steep, possibly arguing against 
the pulsar hypothesis (Fierro et al. 1993), although it should be noted that Halpern et al. (2001) have 
recently presented a strong argument for the association of PSR J2229-I-6114 with 3EG J2227-I-6122, which 
a relatively soft index of 2.24± 0.14. 

Out of 19 sources in Table 1, 12 have the same classification under the two schemes; in Table 2, 6 
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Figure 9: Comparison of the variability indices of the 3EG sources superposed with SNRs and identified 
3EG AGNs. 

out of 7 have the same classification. This confirms, for these sources, that the schemes are statistically 
correlated and that it is safe to consider both indices to smooth out any particular problem with singular 
sources, as apparently is the case for 3EG J0631+0642 and 3EG J0634+0521. 

Figure IHl presents a comparison between the variability indices for those 3EG sources superposed to 
SNRs and the set of 'A' AGNs identified in the same catalog. The distributions are quite different, showing 
a more non-variable population in the case of the sources herein investigated. The most variable sources 
in the set are identified for reference. 

9 Observations and data analysis 

9.1 CO data 

CO is a polar molecule with strong dipole rotational emission in the mm waveband, and is considered a 
reliable tracer of molecular hydrogen, H2, which, though much more abundant, has only a weak quadrupole 
signature. Throughout this paper, molecular gas masses are derived from observations of the J=l-0 rota- 
tional transition of CO at 115 GHz, assuming a proportionality between velocity- integrated CO intensity, 
Wco, and molecular hydrogen column density N(H2). Specifically, we adopt 

N{m)/Wco = 1.8 X 10^° cm-2K-ikm"^s"\ 

the value derived by Dame, Hartmann, & Thaddeus (2001) from an intercomparison of large-scale far- 
infrared, 21 cm, and CO surveys. To estimate total nucleon densities in molecular clouds, we account for 
elements heavier than hydrogen by assuming a mean molecular weight per H2 molecule of 2.76. 

All CO(l-O) data presented here are from the whole-Galaxy survey of Dame, Hartmann, & Thaddeus 
(2001). This survey is a composite of 37 separate CO surveys carried out over the past 20 years with two 
nearly identical 1.2 meter telescopes, one located at the Harvard-Smithsonian Center for Astrophysics in 
Cambridge, Massachusetts, and the other at the Cerro Tololo Interamerican Observatory in Chile. The 
angular resolution of the composite survey is ~ 8.5 arcmin and the velocity resolution ~ 1 km s~^. Many 
of the 37 separate surveys have been published previously and, where appropriate, the earlier papers are 
cited instead of, or in addition to, Dame, Hartmann, & Thaddeus (2001). 

In a few cases, observations of the C0(J=2-1) transition at the same angular resolution, obtained with 
the Tokyo-NRO 60 cm Survey Telescope (Sakamoto et al. 1995), are used to search for enhancements 
of the CO(J=2-1)/CO(J=1-0) ratio as an indicator of SNR-molecular cloud interaction (e.g.. Seta et al. 
1998). 

9.2 Radio continuum data and diffuse background filtering 

In general, the diffuse radio emission of the Galaxy hampers the detection of weak and extended sources 
with low surface brightness, like SNRs. This large-scale diffuse emission can be removed through different 
techniques for astronomical data analysis. These techniques can range from detailed modeling of the 
non-thermal emission in the Galaxy to model-independent filtering algorithms. In this work, many of the 
radio continuum maps have been cleaned of background diffuse contamination using the method originally 
introduced by Sofue & Reich (1979). Basically, the technique consists of convolving the continuum map 
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Figure 10: CO distribution around the remnant IC443 (G189. 1+3.0). The 3EG 7-ray source J0617+2238 
is superposed. Note the positional coincidence of the contours of the latter with part of the most dense 
regions of the CO distribution. The optical boundary of the SNR is superposed as a black contour. The 
optical emission seems to fade in regions where CO emission increases, this indicates that the molecular 
material is likely located on the forground side of the remnant, absorbing the optical radiation. Optical 
contours are from Lasker et al. (1990). 



with a Gaussian filtering beam, producing a new map with a different brightness temperature Tq. A new 
temperature distribution is computed as = T - AT\ for AT^ > and = T for AT^ < 0. In these 
expressions, T is the temperature distribution of the original map, and AT^ = T — Tq are the residuals 
between the original and the convolved maps. The procedure is repeated, now convolving the Tg^ map in 
order to obtain Tq, AT^ and finally Tq. After n iterations, the difference \T'q — Tq~^\ becomes smaller 
than the rms noise, and a map of residuals AT" = T — Tq* is obtained where all diffuse emission with 
size scales larger than the original filtering beam has been removed. The result is completely independent 
of the original mechanism that produced the large-scale emission. In the following sections we apply this 
technique to Effelsberg lOOm-telescope data and to the MOST Galactic plane survey in order to get images 
of the SNRs as clean as possible at low galactic latitudes. 

10 Case by case analysis 

10.1 7-ray source 3EG J0542+2610 SNR G180.0-1.7 

The possible counterparts of this 7-ray source were explored in detail in a recent paper (Romero et al. 2001). 
No known radio pulsar coincides with this source (see Table 3). Additionally, a radio-quiet Geminga-like 
pulsar origin is disfavored a priori because of the high variability and the steep spectral index that this 
source presents (F = 2.67 it 0.22, Hartman et al. 1999). We have searched the EGRET location error box 
for other compact radio sources. Although 29 point-like radio sources were detected, none of them is strong 
enough to be considered a likely counterpart (Romero et al. 2001). The strongest of the sources detected 
have a radio flux one order of magnitude less than those presented by known 7-ray blazars detected by 
EGRET. Moreover, the absence of an X-ray counterpart to this source suggest that it is not an accreting 
source, like a microquasar. 

Some of us suggested that the only object within the 95% error box capable of producing the required 

7-ray flux is the X-ray transient A0535-I-26. This Be/accreting pulsar, not detected at all in the radio 
band, can produce variable hadronic 7-ray emission through the mechanism originally proposed by Cheng 
& Ruderman (1989, 1991). See Romero et al. (2001) for further details. 

On the basis of results discussed in that paper we conclude that 3EG J0542+2610 and G180.0-1.7 
are most likely unrelated. An interesting comparison between this case and one of the EGRET sources 
coincident with the Monoceros Loop is made below. 

10.2 7-ray source 3EG J0617-|-2238 - SNR G189.1-|-3.0 (IC443) 

A detailed description of SNR IC433 was given by Chevalier (1999). We have recently reviewed the 
spatially-resolved multiwavelength spectrum of IC443 and argued that the morphology and spectrum of 
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Figure 11: Hard energy band (3-10 keV) of the IC443 nebula. Most of the thermal emission associated 
with IC443 is not present in this band. The image shows several point sources, besides the plerion nebula 
itself. The nebula can be represented by an ellipse of 8 arcmin x 5 arcmin. The 95% confidence EGRET 
error circle for 3EG J0617+2238 is also shown. From Bocchino & Bykov (2001). 



the 7-ray emission make it a likely hadronic cosmic-ray accelerator (Butt et al. 2002b. 

Seta et al.(1998) have provided an analysis of the CO environment of this remnant. They concluded, 
as was previously reported by Scoville et al. (1977) and Cornett et al. (1977), that 1C433 is interacting 
with several ambient molecular clouds with a total mass of about W^Mq. They also analyzed the ratio 
i?=CO(J=2-l)/(J=l-0) in the environs of IC443 and concluded that parts of the clouds presented an ab- 
normally high value, consistent with shock interaction. The detected value of R exceeds 3 (the average 
galactic value is ~ 0.6) in some regions. Interestingly, the peak of the CO(J=2-1)/(J=1-0) ratio is coin- 
cident with the location of the newly discovered pulsar wind nebula by Olbert et al. (2001), which may 
indicate an alternate way of exciting molecular gas. Dickman et al. (1992) has estimated that the total 
perturbed molecular gas has a mass of 500-2000 Mq. 

In recent years. X-ray observations of IC443 have been carried out. IC443 was a target for X-ray 
observations with HEAO 1 (Petre et al. 1988), Ginga (Wang et al. 1992), ROSAT (Asaoka & Aschenbach 
1994), ASCA (Keohane et al. 1997), and more recently, with Chandra and Beppo-SAX; we discuss the 
latter in more detail below. IC443 was believed to be mostly thermal in the X-ray band (Petre et al. 
1988, Asaoka &: Aschenbach 1994), although it has been discovered to emit hard X-ray emission (Wang 
et al. 1992). Keohane et al. (1997) later found that the hard X-ray emission was localized and non- 
thermal. They concluded that most of the 2-10 keV photons came from an isolated emitting feature 
and from the South East elongated ridge of hard emission. Even more recently, Preite Martinez et al. 
(1999) and Bocchino &: Bykov (2000) reported a hard component detected with the Phoswich Detector 
System (PDS) on BeppoSAX and two compact X-ray sources corresponding to the ASCA sources detected 
with the BeppoSAX Medium-Energy Concentrator Spectrometer (MFCS) (ISAX J0617.1+2221 and ISAX 
J0618.0+2227). ISAX J0617. 1+2221 has also been observed with the Chandra satellite by Olbert et al. 
2001, who also obtained VLA observations at 1.46, 4.86 and 8.46 GHz and a polarization measurement. 
The hard radio spectral index, the amount of polarization, and the overall X-ray and radio morphology led 
them to suggest that the source is a plerion nebula containing a point source whose characteristic cometary 
shape is due to supersonic motion of the neutron star. Bocchino & Bykov (2001) have, in addition, recently 
observed IC443 with XMM-Newton Observatory (see Figure lll|) . They resolve the structure of the nebula 
into a compact core with a hard spectrum of photon index 7 = 1.63^q]^q in the 2-10 keV energy range, 
and found that the nebula also has an extended (~ 8 arcmin x5 arcmin) X-ray halo, much larger than the 
radio emission extension. The photon index softens with distance from the centroid, a behavior also found 
in other X-ray plerions such as 3C58 and G21.5-0.9. Bocchino & Bykov (2001) also looked for periodic 
signals from the NS but found none with 99% confidence level in the 10~^ — 6.5 Hz range. Assuming that 
Lx/E ~ 0.002, the spin-down luminosity of the central object results in E = 1.3 x 10^^ erg s~^. If the 
power-law spectrum of the nebula core region were extrapolated up to the GeV regime, it would provide a 
flux of 2.0 (0.4-15.2) X 10""^ ph cm s , which is consistent with the EGRET flux from the IC443 region 
(see Table 1). However, the fact that the nebula lies outside the 95% confidence circle of the source argues 
against an association. Further, the EGRET source luminosity would require a substantial fraction of the 
estimated spin-down power. 
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Figure 12: Left: CO contours (white) on a background image from the Digital Sky Survey. The map covers 
a somewhat larger region than the left panel Figure, and shows the HII region and young cluster NGC 
2244 producing a hole in the cloud. Most of the Monoceros Loop is also seen faintly. EGRET sources 
contours are marked in black. Right: CO emission plus contours (black) of 1.4 GHz emission to mark the 
Rosette Nebula. The peak of the CO emission if near 3EG J0634+0521. The position of the X-ray source 
SAX J0635+0533 is marked with a star. White lines are the confidence levels of the EGRET sources. 



The H2 mass near SNR IC443 that we report here, 1.1 x lO^M©, is the total molecular mass in the 
IC443 velocity range, v = —40 to +20 km s~^, within a rectangle enclosing the main clump coincident with 
the SNR: /=188.75 to 189.5, b=2.5 to 3.375, see Figure [Till The average density of that region is about 
840 nucleons cm~'^. Assuming that the energy of the explosion was E^i = 0.27, and an unshocked density 
of 0.21 cm~'^ (from a 2D dynamic model by Hnatyk and Petruk 1998), we found that the hadronic fiux 
would be 4.7 x 10~^ photons cm^^ s^^. Indeed, just 10% of the ambient mass is necessary to produce the 
observed flux of 3EG J0617+2238. 

The energy of the explosion and the unshocked density yield a CR enhancement factor, kg ~ 600 
within the SNR, which appears to be unusually high. However, computation of the energy transformed 
into cosmic rays by the direct product A;secR(4/37ri?^) gives 0.4£'5i, a value compatible with that obtained 
for 9 (the efficiency of SN energy conversion to CRs) in the Morfill et al. (1984) prescription discussed 
above. Alternatively, by direct use of Eq. ((T^ . and using the mass considered above and the observed 
7-ray flux, we can obtain an estimate of the enhancement factor kg of 66. The difference between the 
CR enhancement factor of the cloud (/cg = 66) and the SNR (kg = 600) could be explained in a variety 
of ways: most simply, that the CR enhancement predicted by the Morfill et al. (1984) prescription for 
SNRs overestimates the value within the adjacent cloud. Particularly if the cloud abuts the remnant, its 
enhancement will naturally be smaller than that calculated for the SNR interior. As Figure^] shows, the 
optical emission seems to fade in regions where CO emission increases, which perhaps indicates that the 
molecular material is absorbing the optical radiation, abutting the remnant on the near side. The report 
by Cornett et al. (1977) also argues that the molecular mass is located between us and the SNR. In any 
case, we remark that the previous estimations of the CR enhancement factors assume that the explosion 
proceeds in an homogeneous medium, something which we know is not true in this case (Chevalier 1999). 

We note that an electronic bremsstrahlung hypothesis for the origin of the GeV fiux (eg. Bykov et al. 
2000) is difficult to reconcile with the fact that the radio synchrotron emission is concentrated towards the 
rims of the remnant, whereas the GeV source is centrally located (Figure ITU)) . 

It is clear that IC443 will continue to be a primary target for future satellites missions and telescopes. 
A better localization of the EGRET sources, by AGILE, or GLAST, as well as the already approved 
INTEGRAL observations could help much in determining the ultimate nature of 3EG J0617+2238. The 
reader is referred to Butt et al. (2002b) for further analysis of the likely hadronic origin of the 7-ray 
emission. 
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10.3 7-ray source 3EG J0631+0642 and 3EG J0634+0521 - SNR G205. 5+0.5 (Mono- 
ceros nebula) 

The large SNR G205.5+0.5 (Monoceros Loop nebula, 220 arcmin in size) has been thoroughly studied in 
the past. Various papers have proposed that the Monoceros Loop SNR is interacting with the Rosette 
Nebula (e.g., Odegard 1986). A recent study of the stars in NGC 2244, the cluster within the Rosette, finds 
a distance of 1.39ib0.1 kpc (Hensberge et al. 2000); we assume this distance in the following computations. 
In the left panel of Figure \T2\ we show an overlay of the CO contours on an image from the Digital Sky 
Survey (Lasker et al. 1990). It shows very nicely the HII region and young cluster apparently carving a 
hole in the cloud. Most of the Monoceros Loop is also seen very faintly. This justifies the distance adopted, 
under the assumption that the Nebula and the Monoceros Loop are equally distant from Earth. Bloemen 
et al. (1997) presented COMPTEL observations of the Monoceros region, and found excessive 3-7 MeV 
emission which they attributed to nuclear deexcitation lines at 4.44 and 6.13 MeV from accelerated 12 C 
and 16 O nuclei. 

Monoceros was already suggested by Esposito et al. (1996) as a source of 7-rays, and it was also men- 
tioned by Sturner and Dermer (1995) and Sturner, Dermer and Mattox (1996) as a possible case of 7-rays 
production by hadronic interactions. The age of the remnant is not well determined, 3-20 10^ yrs. A study 
of Einstein IPC data (Leahy et al. 1986) shows diffuse X-ray thermal emission in a region corresponding 
to the detection of optical filaments (Odegard 1986). This is only possible if the gas is sufficiently hot, and 
thus if the remnant is sufficiently young. This would be in contradiction with the age one obtains from the 
homogeneous Sedov solutions, which would give an age in excess of 100,000 yrs. One direct interpretation 
of this discrepancy (Leahy et al. 1986) is that the expansion of the remnant proceeds in a non-homogeneous 
multi-component medium, where the homogeneous Sedov solutions are not valid. 

A large region covering the spatial extent of both 3EG sources was studied by Jaffe et al. (1997): 
(198 < I < 214, and —6 < b < 8). They presented an image reconstruction of the region around the 
Rosette Nebula and Monoceros using high-energy (>100 MeV) 7-ray data from EGRET. The resulting 
image showed a 7a extended feature in excess of the expected diffuse emission located at the point-source 
position listed in the EGRET catalog (2EG at that time). These authors proposed that this excess could be 
evidence of an interaction between the Monoceros remnant and the Rossette nebula. They concluded that if 
the 7-ray emission arises solely in the interaction between the two nebulae then the cosmic-ray enhancement 
would be around fes=300. This value appears to be excessively high in this case, should the enhancement 
be the same for all the SNR region. The energy in cosmic-rays, computed using A;secR(4/37ri2^) together 
with the size of the Monoceros remnant (~ 60 pc), imply an energy of the explosion about one order of 
magnitude larger than the assumed E^i ~ 1. This may indicate that the hadronic origin of the 7-rays 
in the interaction of the Monoceros SNR and the Rossete Nebula (i.e. 3EG 3EG J0634+0521) cannot be 
responsible for the entire observed flux. 

Indeed, within the 95% contour of 3EG J0634+0521 there exists an X-ray source SAX J0635+0533, 
and a Be-star/neutron-star X-ray binary pulsar, probably with a relatively short orbital period (Kaaret et 
al. 1999, Cusumano et al. 2000, Nicastro et al. 2000). The hard X-ray source SAX J0635+0533 shows 
pulsations at a period of 33.8 ms but no radio flux was detected at the Be-star position (see below). SAX 
J0635+0533 might be, as in the case of A0535+26 mentioned above, a source of 7-rays through hadronic 
processes. Kaaret et al. (1999) suggested that SAX J0635+0533 and 3EG J0634+0521 are related. One 
fact favoring this physical association is that the SAX satellite has not detected extended emission in the 
region of 3EG J0634+0521, as would be the case if the bulk of the radiation were produced in a SNR 
shock. Additionally, the probability for chance positional coincidence between a Be/X-ray binary and an 
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EGRET source is less than 4% (Kaaret et al. 1999). 

The situation, however, is far from resolved. The EGRET source, for instance, is non-variable as it 
would be expected for a binary with eccentric orbit. Recent results reported by Kaaret et al. (2000), 
comparing observations obtained with BeppoSAX and RXTE separated in time by 2 years, showed that 
the period derivative of the pulsar has a lower bound equal to 3.8 xlO~^'^. This value is 30 times larger 
than values found from accreting neutron stars (Bildsten et al. 1997), and it implies a mass accretion rate 
of 6 xlO~^M0 yr^^ (Kaaret et al. 2000, Bildsten et al. 1997), which far exceeds the expected mass cap- 
ture rate of a neutron star ~ IO^^^A/q yr~^, and is even slightly larger than the average mass loss rate of 
Be-stars (Kaaret et al. 2000). Apparently, this would indicate that the X-ray luminosity does not originate 
in the accretion disk, and argues in favor of SAX J0635-I-0533 being a rotation powered pulsar. In this 
case, the value of P would imply a characteristic age of only 1400 years and a high spin-down luminosity 
of 5 X 10^^ erg s~^, out of which less than 0.05% could make a noticeable contribution to the observed 
7-ray flux (assuming a distance of 4 kpc, Kaaret et al. 1999). Additionally, arguing against an accretion 
origin of the radiation, the derived X-ray luminosity (7.7 x 10^^ (dkpc/4 kpc)^ erg s~^) and magnetic field 
strength (~ 10^ G) are too low in comparison to other Be- X-ray binaries such as A0535+26 (Cusumano 
et al. 2000). 

Very recently, Monoceros was the target of the HEGRA Cerenkov telescopes (see Lucarelli et al. 2001). 
HEGRA observed the Monoceros-Rossete region for about 120 hours, with an energy threshold of 500 
GeV and an angular resolution of 0.1 deg, and mapped a 2 x 2 deg^ region centered in the source SAX 
J0635+0533. The EGRET source 3EG J0634+0521 is also within the field of view. Although the flux 
and spectrum have not yet been officially reported, HEGRA found a tentative excess of counts in four 
different pixels (0.2 x 0.2 deg^) within the 3EG contours; interestingly, none of them coinciding with SAX 
J0635+0533 (Lucarelli et al. 2001). It is possible that TeV emission coming from the binary is being 
re-absorbed in its neighborhood, as in the case studied by Romero et al. (2001). 

What HEGRA observations seem to imply is that the marginally significant TeV radiation has an 
extended origin, different from that producing GeV photons in the binary system SAX J0635+0533, but 
this remains to be confirmed. 

In the right panel of Figure we show the CO emission plus contours of 1.4 GHz continuum emission 
marking the Rosette. Using the standard CO-to-i?2 mass conversion and a mean molecular weight per 
H2 molecular of 2.76, we calculate a total H2 mass for the associated cloud (in the region I = 205 to 
209, b = —3 to —1, and velocity range v = — 5 to 30 km s"^) of 1.2 x 10^ M©. We have calculated the 
molecular masses in small rectangles enclosing the two 3EG sources near the Rosette. Using a distance of 
1.39 kpc we get the following results: In the case of 3EG J0634-I-0521 the region considered is I = 205.5 
to 206.875 and b = -2.125 to -0.375, and the mass is 2.0 xlO^ Mq. For 3EG J0631+0643, in the region 
I = 204.25 to 205.5 and b = -2.125 to -0.75, the mass is 4.7 xlO"^ Mq. In both cases, the velocity range 
considered is f = —5 to 30 km s~^ and the formal error on these masses, based on the instrumental noise, 
is ~ 0.2 X 10^ Mq. 

Because of uncertainties on the nature of this SNR (for instance, the controversy on the SNR age) Mor- 
fill et al.'s method is unreliable for estimating the SNR GeV- flux. However, using directly Eq. (|13|) and 
the observed flux (Table P) , we can estimate the value of fcg needed to generate the observed flux, resulting 
in kg ~ 6.5 for 3EG J0634-I-0521. Because of the high molecular density, just a modest enhancement of the 
cosmic-ray density can explain a substantial part of the detected 7-ray flux. We suggest, then, that 3EG 
J0634-I-0521 might be a composite source: SAX J0635-I-0533 might be responsible for part of the GeV flux. 
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Figure 13: (a) Longitude-velocity map of CO integrated over 1 degree of Galactic latitude roughly centered 
on the SNR G284.3-1.8 (MSH 10-53), b = -2.5 to -1.5. The longitude of the SNR is indicated by the dotted 
vertical line, (b) Spatial map of CO integrated over the velocity range -22 to 3 km s~^. The plus sign 
marks the center position of the SNR G284.3-1.8, whose size is 25 arcmin. The dotted circle is the 95% 
confidence radius about the position of 3EG J1013-5915 (Hartman et al. 1999). Note that the longitude 
range (x axis) of both maps is the same. 



as well as the bulk of the emission at X-ray energies. The interacting SNR and Rossete Nebula might also 
contribute to the flux in the GeV range, and would provide the bulk of the possibly detected TeV emission 
from the region. One direct way to test this scenario would be through an analysis of the spectrum from 
GeV to TeV. In the case of a composite source, there should be a break in the spectrum between GeV and 
TeV energies, the latter corresponding only to the accelerated particles in the SNR remnant. 

In the case of 3EG J063H-0643, a CR enhancement value of just fcg ~ 3 can explain the observed GeV 
flux. New high-sensitivity radio measurements of the region would be of great value in determining the 
relative importance of hadronic and leptonic 7-ray emission. 

10.4 7-ray source 3EG J1013-5915 - SNR G284.3-1.8 (MSH 10-53) 

For this 3EG source, a natural candidate to generate a significant part of the 7-ray emission seems to 
be the recently discovered pulsar PSR J1013-5915; see Table 3 (Camilo et al. 2001). This pulsar has a 
characteristic age of r = 21 kyr and a spin-down luminosity oi E = 2.6 x 10^^ erg s^^. If only the pulsar is 
considered, the efficiency required for converting spin-down luminosity into 7-rays is r/ ~ 0.5% (Camilo et 
al. 2001), well within the range of efficiencies for previously confirmed 7-ray pulsars detected by EGRET. 

Another Parkes' pulsar, PSR J1013— 5934, is also coincident with 3EG 1013—5915, but it can be ruled 
out on the basis of energetic arguments. So also the Princeton pulsar PSR J1012-5857 (Table 3), which 
requires an efficiency of about 100% in the generation of the 7-ray emission from the spin-down losses. It is 
also interesting to note that 3EG J1013— 5915 has a photon spectral index softer than typical for pulsars: 
r = 2.32 ±0.13. 

The remnant is located in the near side of the Carina spiral arm, in a region with a high density of 
molecular clouds. Ruiz & May (1986) found ffiamentary optical emission associated with the remnant. 
These authors also found clear evidence of at least three small CO clouds interacting with G284.3-1.8. 
Other small clouds could have been disrupted by the supernova blast wave and are now forming the shell. 
Their CO(J=1-0) mm line data shows sudden changes with position in radial velocity, and the presence 
of broad asymmetric lines with peak-shoulder profiles, both of which indicate a shock wave disruption of 
the CO clouds. Since Ruiz & May (1986) gave only an upper limit for the mass of the shell, we have 
re-analyzed the gas content for this region. As the longitude- velocity map in Figure shows, nearly all 
CO emission in the general direction of the 3EG source and the SNR G284.3-1.8 lies in the velocity range 
-22 to 3 km s"^ A CO map integrated over this range is shown in Figure 113b . The mean velocity of 
the emission is ~ —9 km s^^, which is consistent with the terminal velocity in this direction, suggesting a 
distance of approximately 2.1 kpc. However, since radial velocity changes very slowly with distance in this 
direction, the uncertainty on the kinematic distance is large, approximately ±1 kpc. We adopt a distance 
of 2.9 kpc, the value inferred by Ruiz & May (1986) based on optical observations of the SNR filaments. 
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Figure 14: CO integrated over the velocity range of the far Carina arm, to 45 km s . Contours: 843 
MHz continuum from the MOST Galactic plane survey (Green 1997); the survey has been smoothed to a 
resolution of 3 arcmin to highlight extended sources. The contour interval is 0.04 Jy/beam, starting at 0.04 
Jy/beam. The dotted circle is the 95% confidence radius about the position of 3EG 1102-6103 (Hartman 
et al. 1999). G290.1-0.8 is the SNR MSH 11-61A (Kirshner & Winkler 1979). Both HII regions are in the 
catalog of Georgelin &: Georgelin (1970) and have recombination line velocities in rough agreement with 
that of the complex. The component clouds A and B are discussed in the text. 



the S — D distance, as well as the CO kinematics. 

There is no CO detected at any velocity toward the nominal center of the 3EG source (/ = 283.93, 
b = —2.34), although the total molecular mass within the 95% confidence radius of the 3EG source (dotted 
circle in Figure IT^ ) is 5.9 x W^Mq. Most of this mass does not coincide with the SNR, which is completely 
included within the 3EG source. Since this emission does not form a single well-defined cloud, it's possible 
that it arises from gas spread over quite a large distance along the line of sight, perhaps 1 — 2 kpc. If so, 
the gas density near the SNR could be quite low. Our study therefore reinforces the idea that it is most 
likely the pulsar, and not the hadronic or bremsstrahlung emission from the SNR neighborhood, that is 
responsible for 3EG J1013-5915. 

10.5 7-ray source 3EG J1102-6103 - SNR G290.1-0.8 (MSH ll-61A)/289.7-0.3 

Sturner & Dermer (1995) proposed that this 7-ray source, in its 2EG J1103-6106 incarnation, may have 
been related to SNR G291. 0-0.1. However, the more precise localization in the 3EG catalog shifted the 
source's position such that it is no longer superposed with that SNR. 

Zhang & Cheng (1998) argued against a newly discovered young radio pulsar, PSR J1105-6107 (Kaspi 
et al. 1997), as the source of the observed high-energy 7-ray emission. Its age (~ 6.3 x 10^ yr) also seems 
high for a Vela-like pulsar. In addition, the photon spectral index is very soft, F = 2.47 it 0.21, though 
this in itself does not disqualify a possible pulsar origin, as seen in the case of PSR J2229+6114/3EG 
J2227+6122 (Halpern et al. 2001b). 

The line of sight to 3EG 1102-6103 intersects both the near and far sides of the Carina spiral arm, at 
velocities near —20 km s~^ and +20 km s^^ respectively. There is a distinct gap in the near side of the 
Carina arm in the direction of the 3EG source, with almost no molecular gas within ~ 1 deg of the source 
direction (see, e.g.. Figure 2 of Dame, Hartmann, & Thaddeus 2001). On the other hand, as Figure ITU 
shows, there is a very massive molecular complex in the far Carina Arm overlapping the direction of the 
3EG source; this complex is No. 13 in the Carina Arm cloud catalog of Grabelsky et al. (1988). There is 
little doubt that the two component clouds labeled A and B in Figure are part of the same complex, 
since they have approximately the same velocity of 22 km s~^, and are connected smoothly by weaker 
emission, also at the same velocity. Also, the HII regions are evidence of abundant on-going star formation 
in this molecular complex which additionally supports the association of the SNR. Assuming a flat rotation 
curve beyond the solar circle, the kinematic distance of the complex is 8.0 kpc and its total molecular mass 
is 2.1 xW^Mq. 

It is worth noting that the composite CO line profile of cloud B is very broad and complex, suggesting 
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possible interaction with SNR G290. 1-0.8. In the case of cloud A, its radius (~48 pc) and composite 
linewidth (17 km s~^ FWHM) are roughly consistent with the radius-linewidth relation found for large 
molecular complexes by Dame et al. (1986). For cloud B, however, its linewidth (~27 km s~^) is about a 
factor of 3 too large compared to its radius (~28 pc). We can also see that the coinciding SNR G289. 7-0.3 
is far from Cloud B, in a region of low molecular density. It is extremely unlikely that this SNR is related 
with the 3EG source in question. The only remaining candidate is, then G290. 1-0.8. The total molecular 
mass within the 95% confidence radius of the 3EG source (dotted circle in Figure ITU) is 7.7 xIO^Mq and 
most of it is localized in Cloud B (4.5 xlO^M©). 

Assuming typical values for the energy of the explosion (£51=1) and the unshocked ambient density 
(n = O.lcm"'^) we obtain a CR enhancement factor of ~ 250. Assuming that the same CR enhancement 
is applicable to the cloud overpredicts the EGRET flux by about a factor of 10. Thus, it is likely that 
the average CR enhancement factor within the cloud is ten times lower than within the SNR, a reason- 
able result. It is possible, then, that 3EG J1102-6103 and SNR G290.1-0.8 are indeed related. Note 
that Bremsstrahlung, which we have neglected here, will contribute still more to the predicted flux from 
SNR-cloud interactions. If the outlined scenario is correct, GLAST and AGILE ought to observe a strong, 
compact 7-ray source coincident with the position of Cloud B. 

An alternative, and promising, hypothesis for explaining the high-energy emission is stellar winds 
collisions, as developed by Eichler & Usov (1993) and Benaglia & Romero (2003). Recently, Contreras 
et al. (1997) have provided convincing evidence for non-thermal radio emission from the colliding winds 
region in the stellar system Cygnus 0B2 No 5.^ 

The position of their radio-imaged shocked region is consistent with the inferred location of the contact 
discontinuity of the wind-wind interaction of the constituent stars. Benaglia et al. (2001) have argued that 
the source 3EG J2033+4118 could be mainly due to inverse Compton scattering of the stellar photons by 
the locally accelerated electrons. Similarly, in the present case, 3EG J1102-6103 might be the result of 
7-ray production by the stellar winds of the early-type stars WR37, WR38, WR38B and WR39, all located 
within the 95% confidence contour of the 7-ray source and at 2 kpc from the Sun. In particular, WR39 
presents an unusually strong wind with a terminal velocity of about 3600 km (Romero et al. 1999a and 
references therein). Non-thermal radio emission at the mJy level has been recently detected at ~3 arcsec 
from the optical position of the star by Chapman et al. (1999). This emission is a clear indication of the 
existence of a population of relativistic electrons in the region. Chapman et al. (1999) have suggested that 
particle acceleration could be occurring at the region where the wind of WR39 collides with the wind of the 
neighboring Wolf-Rayet star WR38B. This hypothesis is supported by the fact that synchrotron radiation 
is located between both stars (2 arcsec from WR38B). The relativistic electrons should interact with UV 
photons from the star, producing IC 7-rays that could explain part of the emission of 3EG J1102-6103. 

Fortunately, the peak of the spectral energy distribution should be in the IBIS energy range, the Im- 
ager on-Board INTEGRAL. Using Benaglia et al.'s (2001) model with a spectral index of —2, we get an 
integrated fiux for the energy interval 100 keV - 200 keV of 1.2 x 10"^ ph s'^ cm . For the entire IBIS 
energy range (20 keV-10 MeV) the expected value is 1.2 x 10^^ ph s~i cm . In addition to this wind-wind 
contribution, single stars might also be sources of 7-rays in the IBIS energy band through the IC emission 

''The HEGRA Cherenkov telescope array group recently reported a steady and extended unidentified TeV gamma-ray 
source lying at the outskirts of Cygnus OB2, the most massive stellar association known in the Galaxy, estimated to contain 
2600 OB type members alone (Aharonian et al. 2002). Butt et al. (2003) reported on near-simultaneous follow-up observations 
of the extended TeV source region with the CHANDRA X-ray Observatory and the Very Large Array (VLA) radio telescope. 
The broadband spectrum of the TeV source region favors a predominantly nucleonic rather than electronic origin of the 
high-energy flux, possibly in a way similar to that proposed by Romero and Torres (2003) in the case of NGC 253. 
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Figure 15: CO integrated over the velocity range of the Centaurus arm, -75 to -35 km s . Contours: 843 
MHz continuum from the MOST Galactic plane survey (Green et al. 1999); the survey has been smoothed 
to a resolution of 3' to highlight extended sources. The contour interval is 0.01 Jy/beam, starting at 0.01 
Jy/beam. The dotted circle is the 95% confidence radius about the position of 3EG J1410-6147 (Hartman 
et al. 1999). Many of the other radio sources in the map are unrelated HII regions discussed by Caswell 
& Barnes (1985). 



of electrons locally accelerated in shocks at the base of the winds. These shocks are produced by line-driven 
instabilities (see Chen &: White 1991 for a discussion of 7-ray emission from single stars, and Benaglia et 
al. 2001 for a particular application). 

Consequently, further study of the possible association of 3EG J1102-6103 with the SNR G290. 1-0.8 
(MSH 11-61A) is of the utmost importance, since there are at least two scenarios (aside from the pulsar 
possibility) that might well contribute to the observed 7-ray flux. 

10.6 7-ray source 3EG J1410-6147 - SNR G312. 4-0.4 

Although there is no correlated Princeton pulsar within the contours of 3EG J1410-6147, there are two 
Parkes pulsars near the 7-ray source (see Table 6 above). Both of them require an unreasonable efficiency, 
at their dispersion-measure distances, to explain the observed 7-ray flux. Both pulsars are located (at 
least in projection) well within the boundaries of the incomplete shell of SNR G312.4— 0.4 (Caswell &: 
Barnes 1985), to which Yadigaroglu & Romani (1997) estimate aT, — D distance of 1.9 kpc, whereas Case 
&: Bhattacharya (1999) flnd 3.1ibl.O kpc. At either of these distances, the required efficiencies would be 
substantially smaller. The photon spectral index of the 7-ray source, T = 2.12 it 0.14, seems to be in 
the range of other pulsar cases. More recently, Doherty et al. (2002) have provided important new HI 
absorption measurements toward SNR G312.4-0.4 which indicate that it may be much further away, 8.1 
kpc (see below). 

Case & Bhattacharya (1999) have made an in depth study of the possible association between the rem- 
nant and the 7-ray source and concluded that the 7-ray data alone cannot at present provide conclusive 
evidence to decide whether the 7-ray emission from 2EG J1412-6211 is due to a pulsar or SNR-molecular 
cloud interaction, or both. They suggested that CO observations of the environment surrounding G312.4- 
0.4 would help in determining whether a molecular cloud of sufficient mass is present in the right location 
to produce the observed 7-ray intensity. Such observations are presented here. 

The CO emission toward 3EG 1410-6147 is extremely bright and complex, arising mainly from the 
tangent region of the Centaurus spiral arm at < — 30 km s~^. The line of sight also intersects the near 
side of the Carina arm at less negative velocities and the far side of Carina at positive velocities. The cloud 
in the general direction of the 3EG source is also by far the most massive and dense. This cloud, with a 
velocity of —49 km s~^, is labeled A in Figure El The near and far kinematic distances of the cloud are 
3.3 kpc and 8.1 kpc respectively (Clemens 1985); we will adopt the far kinematic distance since it agrees 
with the HI absorption measurements of Doherty et al. (2002) towards SNR G312.4-0.4. 

The mass of Cloud A is quite uncertain owing to the fact that molecular gas in both the near and far 
sides of the Centaurus arm probably contribute to the CO emission at the cloud velocity; the high-velocity 
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Figure 16: Relative positions of 3EG J1639-4702 (which occupies all the box) and the SNRs G333. 8-0.1, 
G338. 1+0.4, and G338. 3+0.0, which are very small in comparison. The superposed grey-scaled levels show 
the radio emission at 843 MHz of the three SNRs as reported in the MOST Catalog prepared by Whiteoak 
and Green (1996). 



Figure 17: CO integrated over the velocity range -65 to -45 km s~^. Contours: 843 MHz continuum from 
the MOST Galactic plane survey (Green et al. 1999); the contour interval is 0.2 Jy/beam, starting at 0.1 
Jy/beam. The dotted circle is the 95% confidence radius about the position of 3EG J1639-4702 (Hartman 
et al. 1999). 



limit of the cloud is also uncertain owing to blending with emission at higher (less negative) velocities. If 
we assume that all of the emission in the velocity range of Figure IHK- 75 to -35 kms"^) arises from cloud 
A at 8.1 kpc, the total molecular mass within the 95% confidence radius of the 3EG source (dotted circle 
in Figure 15) is 3.3 xIO^Mq; given the uncertainties just discussed, the actual mass might be as much as 
a factor of 2 lower. 

Even at such a distance, the large quantity of molecular material is sufficient to explain a significant 
part of the 7-ray flux observed. As in the case of 3EG 1102-6103, with the usual assumptions, a CR 
enhancement factor of ~100 is necessary to explain the bulk of the 7-ray emission from 3EG J1410-6147. 
It is most likely, however, that this EGRET source is a composite. New X-ray observations could be used 
to study the pulsars and evaluate their 7-ray emissivities. An extrapolation of the GeV spectrum to the 
TeV regime, if there is no break, would give a flux of 1 xl0~^'^ erg cm~^ s~^, which is above the HESS 
sensitivity in the range 500 GeV-10 TeV. Even with a substantial break in the spectrum the 3EG source 
should be observable by HESS (see the study on SNR TeV observability below. Section n3|) . 

10.7 7-ray source 3EG J1639-4702 - SNR G337.8-0.1/338.1+0.4/338.3+0.0 

Figure El shows the relative positions of these SNRs within the large location contours of 3EG J1639- 
4702. One Princeton pulsar is within the contour of 3EG 1639-4702, but it can be ruled out as a possible 
counterpart because of the required energetics. In addition, three Parkes pulsars coincide with the same 
3EG source (see Table 3). Two of them can be immediately discarded based on the same grounds: the 
required efficiencies are unphysically high. However, PSR J1637-4642 seems to be a plausible candidate. 
Only a 12% efficiency would be needed to convert this pulsar into a plausible counterpart for the origin 
of the 7-ray emission. Although the spectral index, F = 2.50 ± 0.18, seems quite soft in comparison with 
detected EGRET pulsars, the work of Halpern et al. (2001) suggests that a soft spectral index does not 
automatically rule out a pulsar origin of the 7-rays: they present a strong case for PSR J2229+6114 being 
responsible for 3EG J2227+6122, even though it has a high-energy spectral index of 2.24+ 0.14. 

Based on HI absorption seen all the way up to the terminal velocity, Caswell et al. (1975) placed the 
SNR G337. 8-0.1 beyond the tangent point at 7.9 kpc. Koralesky et al. (1998) detected maser emission in 
the SNR at —45 km s~^, implying a far kinematic distance of 12.4 kpc. As Figure [T71 shows, there is a 
very massive giant molecular cloud adjacent to the SNR in direction and close to the associated maser in 
velocity (—56 km s~^). The far kinematic distance for this giant molecular cloud is favored by (1) its likely 
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Figure 18: Intensity map of the CO(J=1-0) transitions in the region around RX J1713. 7-3946, from Butt 
et al. (2002a). Two massive clouds, called Cloud A and Cloud B, are indicated. The X-ray contours of the 
SNR (Slane et al. 1999) are superposed in black, as well as the location confidence contours of the GeV 
7-ray source 3EG J1714-3857 (coincident with Cloud A) and the significance contours of the TeV detection 
of the remnant by Enomoto et al. (2002), mostly coincident with the X-ray radiation. 



association with both the far-side maser just mentioned and a group of far-side HII regions (group 5 in 
Georgelin & Georgelin 1976); (2) its location very close to the Galactic plane; and (3) the radius linewidth 
relation for giant molecular clouds (Dame et al. 1986). The mean velocity of the complex is —56 km s^^, 
implying a far kinematic distance of 11.8 kpc. 

The giant molecular cloud has recently been discussed by Corbel et al. (1999) because an adjacent 
cloud (just outside the velocity integration range of Figure ITT)) apparently harbors the soft 7-ray repeater 
SGR 1627-41. Corbel et al. suggest that collision or tidal interaction between these two giant molecular 
clouds may have set off the burst of star formation evident in both. 

Taking the total CO luminosity of the giant molecular cloud to be that in the range I = 337.625 to 
338.25, b = —0.25 to 0.25, and v = —65 to —45 km s~^, the total molecular mass is 5 x 10^ M©; this mass 
may be overestimated by 10-20% owing to the inclusion of emission from gas at the same velocity at the 
near kinematic distance. Even with this correction, this giant molecular cloud ranks among the few most 
massive GMCs in the Galaxy (see, e.g.. Dame et al. 1986); its composite CO linewidth of ~ 20 km s~^ is 
correspondingly very large. Adopting a mean radius of 0.31 deg, or 65 pc at 11.8 kpc, the mean nucleon 
density of the cloud is 176 cm~'^. The total mass within the 95% confidence radius of the 3EG source is 
7.6 X IO^Mq; this mass too may be overestimated by 10-20% owing to inclusion of near-side emission. 

The enhancement factor obtained from Eq. I|15() is large for typical parameter values. However, since 
the SNRs seems to be immersed in the molecular cloud, the Morfill et al. (1984) prescription may be an 
oversimplification in this case. This case is similar to that of 3EG J1903+0550 in that we have a distant 
SNR that we would not expect to be able to detect with EGRET in the neighborhood of a very large 
molecular cloud. AGILE observations, in advance of GLAST, would greatly elucidate the origin for this 
3EG source, since even a factor of 2 improvement in resolution would be enough to favor or reject the SNR 
connection. 

10.8 7-ray source 3EG J1714-3857 - SNR G348.5+0.0/348.5+0.1/347.3-0.5 

The supernova remnant RX J1713. 7-3946 is probably the most convincing case for a hadronic cosmic-ray 
accelerator detected so far in the Galaxy. Butt et al. (2001) noted the positional coincidence of the nearby 
7-ray source 3EG J1714-3857 with a very massive ('^ 3 x 10^ MQ) and dense (~ 500 nucleons cm~^) 
molecular cloud that is clearly interacting with the SNR RX J1713.7-3946 (G347.3-0.5) (Slane et al., 1999; 
Butt et al., 2001). Figure IT51 shows the CO(J=1-0) line intensity distribution in the vicinity of the SNR. 
The remnant is a strong X-ray source (Slane et al. 1999) whose ROSAT contours are indicated in the 
figure. Two massive clouds, which we called Cloud A and Cloud B, can be seen. The first one is coincident 
with the 7-ray source 3EG J1714-3857, whose location confidence contours are also superposed in the 
figure. The X-ray emission is produced by TeV-range electrons radiating by the synchrotron mechanism in 
the local magnetic field. These same electrons were suggested to be responsible, through IC up-scattering 
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Figure 19: The radio synchrotron spectrum which would be expected from the region of the shocked molec- 
ular material located towards the NE of the remnant RX J1713. 7-3946 (Cloud A), under the assumption 
that the observed GeV flux were due to electron/positron Bremsstrahlung. Since this spectrum violates 
the upper limit (dark) derived from the non-detection of the cloud in the radio band by a factor of ~ 20 
at 843 MHz (Slane et al., 1999), we can rule out a predominantly leptonic origin of the GeV luminosity. 
Furthermore, if the GeV flux of 3EG J1714-3857 were of electronic origin, the cloud region would outshine 
even the radio-brightest NW rim of the remnant, which is found to be emitting only 4ibl Jy at 1.36 GHz 
(Ellison et al. 2001), as shown by the light data point. An assumed low frequency turnover at ~100 MHz 
in the radio spectrum is shown by the red dotted line. 



of cosmic microwave background photons, for the TeV 7-ray emission detected from the NW rim of the 
remnant by Muraishi et al. (2000) and Butt et al (2002a). 

We have previously provided measures of line intensity ratios R =CO(J=2-1)/(J=1-0) for the entire 
region demonstrating that the SNR is likely interacting with Cloud A (Butt et al. 2001). For this cloud, 
R ~ 2.4, more than 3.5(T above the average Galactic value. 

Upper limits to the continuum radio emission of Cloud A and the use of Eq. H25|) above, allow us to rule 
out a Bremsstrahlung origin of the GeV radiation from 3EG J1714-3857 (see Figure [T9|) : the electron flux 
needed to explain the GeV source in terms of bremsstrahlung emission overpredicts the radio synchrotron 
emission for any reasonable molecular cloud magnetic fields (Cructher, 1988,1994,1999). The GeV 7-rays 
seem to be the result of 7r*^-decays produced when the population of cosmic-rays accelerated at the remnant 
shock are injected into the dense medium of Cloud A. We estimate a cosmic-ray enhancement factor in the 
range 24 < /cg < 36 given the parameters of the SNR (Slane et al. 1999). 

In addition, the 7-ray spectrum of 3EG J1714-3857 (Hartman et al. 1999) is consistent with a narrow 
spectral bump at ~ 70 MeV that could correspond to the signature of the pion-decays resulting from an 
enhanced population of low energy (E~l GeV) protons. This apparent peak, although highly suggestive, 
is not statistically significant and improved observations are needed to confirm its existence. 

Uchiyama et al. (2002a) have reported the discovery of extended (10' x 15') and hard (spectral shape 
described by a flat power-law photon index F = I.0IQ 3) X-ray emission from the position of Cloud A, using 
ASCA data. This emission is interpreted as bremsstrahlung from a Coulomb-loss-flattened distribution of 
nonthermal low-energy protons in the cloud or mildly relativistic electrons (see also Uchiyama et al. 2002b) 
. Uchiyama et al. (2002a) estimate that the energy content in subrelativistic protons within the cloud far 
exceeds that in the relativistic protons, say by a factor ~ 80. The explanation could be that the bulk of 
the more energetic particles have already diffused from the cloud whereas the sub-relativistic population 
is captured there. Alternatively, energetic secondary leptons may also be producing low-level non-thermal 
X-ray and radio emission in the clouds. 

Regarding the highest energy particles produced in the SNR, a CANGAROO re-observation of the 
NW-rim of RX J1713.7-3946 with a new 10-m reflector (CANGAROO II) has allowed a determination of 
the TeV 7-ray spectrum, which can be fitted with a power-law of photon index F ~ —2.8 (Enomoto et al. 
2002). Such a steep spectrum is hard to explain by IC emission and Enomoto et al. (2002) have claimed 
that the TeV emission is also of hadronic origin. If this were the case, however, the GeV gamma-ray flux 
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Figure 20: Upper panel: MOST image of the SNR G355. 6+0.0 at 0.843 GHz (Gray 1994). The grey-scale 
representation ranges from 0.4 x 10~^ to 10 x 10~^ Jy beam^^. Radio contours are shown in steps of 
1 Jy beam~^, the resolution is 43 arcsec. Part of the 7-ray probability contours of 3EG J1734-3232 are 
superposed. Lower panel: Detailed radio image of the SNR at the same frequency. Radio contours are 
shown in steps of 1 x 10"^ Jy beam~^, starting from 0.5 x 10^^ Jy beam~^. 



would be much higher than observed, as recently noted independently by Reimer & Pohl (2002) and by 
Butt et al. (2002a). (It is conceivable that protons with a hard index of ~-1.8 may alleviate the discrepancy 
but then it is more difficult to explain how such a population of protons could produce the measured -2.8 
index TeV 7-rays) ^). 

The SNR RX J1713. 7-3946 is perhaps the best natural laboratory available today for studying the 
acceleration and diffusion of cosmic-rays. The unique combination of a relatively close SNR and a group 
of well defined molecular clouds in its surroundings, none of them in front of the remnant itself, makes this 
source a priority target for the forthcoming generation of high-energy instruments such as HESS, AGILE, 
INTEGRAL, and, GLAST, as well as for infrared, radio, mm and sub-mm observatories. 

10.9 7-ray source 3EG J1734-3232 - SNR G355. 6+0.0 

The shell type supernova remnant G355. 6+0.0 (FigureOO]) was first identified in the MOST Galactic Center 
survey (Gray 1994) as a compact ('^O.l deg x ~0.1 deg) radio source with both thermal and non-thermal 
emission. The western limb shows possible indications of an interaction with the ambient diffuse thermal 
gas present at that location. X-ray emission from this SNR has also recently been reported by the ASCA 
X-ray satellite under the designation AX J173518-3237 (Sugizaki et al., 2001). 

The physical relation of even part of the 7-ray flux of 3EG J1734-3232 with SNR G355. 6+0.0, however, 
is unclear because of the lack of information about both the SNR and its environs. The lack of a 7-ray 
spectral index for 3EG J1734-3232 (Hartman et al. 1999, see Table 1) further complicates any attempt to 
connect the SNR with the 7-ray emission. The 7-ray error box is also coincident with a very young open 
cluster, NGC 6383, {l,b) = (355.66,0.05), which is centered around the bright spectroscopic binary HD 
159176 (07V+07V) (eg. van den Ancker et al., 2000). Together with NGC 6530 and NGC 6531, NGC 
6383 belongs to the Sgr OBI association. The nearby radio source G355.3+0.1, also within the 7-ray error 
box, is most likely an HII region at a distance of ~10 kpc (Crovisier et al. 1973). 

Although the Third EGRET catalog lists GeV 1732-3130 (Lamb & Macomb 1997) as an alternate name 
for this source, the large positional offsets indicate that these two may be separate 7-ray sources (Roberts 
et al. 2001a). Interestingly, "bridging" these two sources is a "possibly variable" COS-B source, 2CG 
356+00, located at {l,b) = (356.5, +0.3) (Swanenburg et al. 1981), which may be related to one or both 
of them. 

The report of a bright, transient hard X-ray source, KS/GRS 1730-312 {l,b = 356.6, +1.06) (Vargas 
et al. 1996), within the 95% error ellipse of GeV 1732-3130 may also explain part of the detected 7-ray 
emission from this region. It is possible that in the quiet state KS/GRS 1730-312 was actually seen in 

■'One possible solution would be to introduce the effects of difussion, if the 7-rays originate in one of the nearby clouds and 
not in the SNR itself as suggested by CANGAROO team. See also Uchiyama et al. 2002c. 
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ASCA data, as "src 1" in Roberts et al. (2001). A mild indication of variability for 3EG J1734-3232 
(I = 2.9, r = O.OOg 

!oo) appears within the /-scheme, and would support the flaring hard X-ray source being 
associated with it, although this is not conclusive with the data now at hand. The transient source could 
also be separately connected to GeV 1732-3130/2CG 356+00. 2CG 356+00 is listed as "possibly variable", 
supporting such an argument (Swanenburg et al. 1981). 

It is clear that no conclusive determination can be made regarding possible association among SNR 
G355.6+0.0, GeV 1732-3130, 2CG 356+00, and 3EG J1734-3232 until the sizes of the 7-ray error boxes 
are significantly reduced by future GeV telescopes. It is almost certain that many GeV sources, especially 
those located towards the inner Galaxy, will eventually be resolved into several separate sources. 

10.10 Near the Galactic Center: 7-ray source 3EG J1744-3011 - SNR G359.0- 
0.9/359.1-0.5 and 7-ray source 3EG J1746-2851 - SNR GO.0+0.0/0.3+0.0 

Since the 3EG sources J1744-3011 and J1746-2851 lie near the very confused Galactic Center region, an 
analysis of the sort we present for other SNR-EGRET source pairs is not possible in this case. Galactic 

Center region sources must to be considered as possibly extended (most likely composite) and confused, 
embedded in high and structured background, and the analysis procedures used for other sources in the 
EGRET catalogs may not apply (Mayer-Hasselwander et al. 1998). 

A detailed discussion of the Galactic Center region is beyond the scope of this report; for this the reader 
is referred to Markoff et al. (1999), Yusef-Zadeh et al. (2000), Meha & Falcke (2001), as well as the book 
edited by Falke et al.(1999). Very recently, Goldwurm (2001) has presented a review of the high-energy 
emission detected from the direction of the Galactic center; the reader is referred to that paper for details 
on different X-rays observations of the Galactic Center region. 

A complete discussion of EGRET observations of the Galactic Center region was presented by Mayer- 
Hasselwander et al. (1998). They found 5488+516 counts above 30 MeV, representing a high significance 
excess. The region analyzed was larger than the position of the central EGRET sources by several degrees 
(see Figure 1 in Mayer-Hasselwander et al. 1998). They also re-analyzed COS-B data, showing that the 
COS-B observations of the region were in agreement, despite previous claims, with the more recent EGRET 
data. Several objects in the region are potential counterparts for the 7-ray radiation detected, including, 
for example, GRO 1744-28, 2S1743-2941, E1740-2942, PSR 1742-30, GRS 1736-297, GRS 1739-278, and 
GX 359+02. However, not all of these coincide with the positions of the 3EG sources we are considering. 
The 7-ray fluxes in different energy bands give only the hint of variability, which is consistent with 3EG 
Catalog estimates for / and r (Table 1). The photon spectrum of the diffuse emission of the Galactic 
Center region shows a clear break at energies about 1 GeV, with a significant steepening thereafter (it 
shifts from —1.3 to —3.1). The hard spectrum at energies above 100 MeV has to be compared with the 
already hard values obtained with standard EGRET analysis techniques, quoted in Table 1 for the sources 
of interest. These hard values argue against 7-rays being produced in diffusive processes, i.e., by ambient 
matter-cosmic-ray interactions. However, this is yet to be confirmed. 

Yusef-Zadeh et al. (2002) have suggested that the central 7-ray source may be due to the interaction 
of the GO. 13-0. 13 molecular cloud with the diffuse and filamentary X-ray features discovered using Chan- 
dra, all lying within the 95% confidence location contours of 3EG J1746-2851. The hard spectrum of the 
EGRET source (—1.7) seems to match the cloud spectrum at about 10 keV when extended down to X-ray 
energies. Electron Bremsstrahlung and inverse Compton may be responsible for the GeV emission. Pohl 
(1997) raised the possibility that the radio arc at the Galactic Center could be the counterpart of the 
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Figure 21: Relative positions of 3EG J1800-2338 and SNR W28. Contours in step of 5 mK, starting from 
15 mK, from the 2695-MHz map obtained with the Effelsberg 100-m single dish telescope (Fiirst et al. 
1990). 



high-energy 7-ray source. Existing radio data on the arc support the view that its synchrotron emission 
originates from cooling, initially monoenergetic electrons that diffuse and convect from their sources to 
the outer extensions of the arc. If the source of high-energy electrons coincides with the Sickle region 
(GO. 18-0.04), as indicated by the radio data, then the ambient far-infrared photons could be subject to 
inverse Compton interaction by high-energy electrons. Pohl (1997) showed that the predicted 7-ray emis- 
sion depends mainly on the magnetic field strength in the arc and that both the flux and the spectrum of 
the central source could be explained by such a process. 

On the other hand, the starving state of accretion flow around the supermassive black hole make it a 
dubious counterpart for the high-energy radiation. Note, however, that from an statistical point of view 
the probability of such a good agreement in the positions of the Galactic Center and the 3EG source is 
about 10~^ (Mayer-Hasselwander et al. 1998). However, early scenarios were presented in which the 7-ray 
emission is produced by the wind accretion from the nearby IRS16 cluster (Melia 1992, Mastichiadis & 
Ozernoy 1994). More recently, Markoff et al. (1997) were able to reproduce the observed spectrum with a 
combination of synchrotron radiation and pion decay. If the 7-ray flux is directly related to the dissipation 
of gravitational energy, i.e. if it is produced by relativistic particles energized by a shock within the infalling 
plasma, Sgr A* could still be the source of the 7-rays observed. However, in a refined analysis, Markoff et 
al. (1999), using data from the 3EG catalog and an improved physical treatment, concluded that this was 
not the case. 

Forthcoming satellites, particularly INTEGRAL, will scrutinize the Galactic center region and hopefully 
resolve the nature of this interesting region. 

10.11 7-ray source 3EG J1800-2338 - SNR G6.4-0.1 (W28) 

This EGRET-SNR positional coincidence was originally proposed by Sturner et al. (1996) and Esposito et 
al. (1996). However, the previous 2EG J1801-2321 source does not coincide with 3EG J1800-2338, since 
the new position has shifted position by about half a degree. W28 was also presented as a possible candi- 
date for an association with a COS B source by Pollock (1985), after noticing that SNR W28 appears to 
be interacting with molecular clouds. Dubner et al. (2000) have made a recent study of the remnant using 
the VLA. They concluded that the remnant is indeed interacting with molecular clouds in the vicinity, 
and observed maser emission, as earlier reported by Claussen et al. (1997, 1999). Arikawa et al. (1999) 
observed W28 and mapped the CO(J=3-2) and CO(J=1-0) rotational transition lines toward the remnant, 
and also concluded that the remnant is interacting with the clouds. The mass of the clouds was found to 
be 2 X 10^ M0. This mass, however, is with respect to the 2EG source position. A new molecular mass 
estimate with respect to 3EG J 1800-2338 is given below. 

The H-alpha filaments in W28 have a mean velocity of 18ib5 km s^^ (Lozinskaya 1974). HI absorption 
measurements by Radhakrishman et al. (1972) are consistent with this velocity, since absorption features 
at 7.3 and 17.6 km s~^ are seen against the SNR continuum. As Figure l22l (lower panel) shows, there is 
a molecular cloud at about the same velocity (~19 km s~^) which is most likely the birth place of the 
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Figure 22: (a) CO integrated over the velocity range to 28 km corressponding to the distance of the 
SNR W28. Contours: 4850 MHz continuum from the survey of Condon et al. (1991); the lowest contour is 
at 0.32 Jy/beam and the contours are logarithmically spaced by a factor 2. The three labeled sources are 
HII regions: M20 the Trifid Nebula and M8 the Lagoon Nebula. The dotted circle is the 95% confidence 
radius about the position of 3EG J1800-2338 (Hartman et al. 1999). (b) CO intensity integrated over the 
latitude range of the molecular gas associated with W28, b = -0.5 to -0.125 deg. 



supernova progenitor. Directly toward the radio-bright rim of W28 ( I ~6.6 deg, b ~ -0.3 deg), which 
has been proposed by Wootten (1981) and Arikawa et al. (1999) as the site of SNR-giant molecular cloud 
interaction, a jet-like CO feature is seen in Figure 122 extending to a smaller cloud at ~7 km s~^. Arikawa 
et al. (1999) proposed the 7 km component as the systemic velocity of W28, whereas Claussen et al. 
(1997) suggested that the entire 7 km s^^ cloud has been accelerated by the SNR from the main cloud 
~10 km s~^ higher in velocity. The alignment of the 7 km s~^ cloud with the bright interacting rim of 
W28 and the jet-like feature linking it to the larger cloud at higher velocity supports the Claussen et al. 
proposal. 

The low longitude of W28 (6.5 deg) makes both the kinematic distance and the mass of the associated 
molecular cloud difficult to measure. Assuming a systemic velocity of 19 km s~^, the near kinematic dis- 
tance is 3.7ibl.5 kpc. If all emission in the velocity range — 28 km s~^ is associated with W28 (see Figure 
122b ). the total molecular mass within the 95% confidence radius of the 3EG source is 3.9 x 10^ Mq. Owing 
to the severe velocity crowding at this Galactic longitude, this mass estimate may be overestimated by as 
much as 50%. 

Even reduced by 50%, the ambient molecular material can still account for the observed EGRET flux. 
Ignoring bremsstrahlung, a CR enhancement factor of about ks ~ 20, similar to that found in other cases, 
is necessary to explain the emission by hadronic interactions. This result is stable against reasonable vari- 
ations in the input parameters and is compatible with consistency tests. Note that the distance used here 
is near the upper limit of those shown in Table [IJ a smaller distance would make a physical association 
even more likely. Velazquez et al. (2002), in a recent analysis of large-scale neutral hydrogen around W28, 
adopt a distance of ~1.9 kpc. They concluded that the SN energy was ~ 1.6 X 10^° erg about 3.3 x lO'' 
years ago. 

An intriguing possibility in this case is that there are actually two (or more) 7-ray sources, each associ- 
ated with a different molecular cloud in Figure and /or with the pulsar PSR B1758-23 (which, coincides 
with the GeV source as reported in the Roberts et al. (2001) catalog). Investigation of this possibility 
will require the next generation of GeV and TeV-telescopes, with their improved sensitivity and angular 
resolution. 

W28 is one of the few remnants observed with Cerenkov telescopes. Rowell et al. (2000) observed 
it with CANGAROO 3.8-m telescope and were able to set an upper limit on the flux (for photons with 
E > 1.5 TeV) of a diffuse source encompassing the clouds discovered by Arikawa et al. (1999) and part of 
the 3EG source: 6.64 xlO~^^ ph cm~^ s~^. A simple extrapolation of the 3EG flux, with the same spectral 
index, up to TeV energies yields a value higher than this upper limit by more than an order of magnitude 
(see Figure 5 of Rowell et al. 2000). This implies the existence of a break in the spectrum in the GeV- TeV 
region. But even considering such a break (see Section 1131 and Table EB the 3EG J1800-2338 region could 
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Figure 23: High-resolution radio map of the nearby star LS 5039 obtained with the VLBA and the VLA 
in phased array mode at 6 cm wavelength. The presence of radio jets in this high- mass x-ray binary is 
the main evidence supporting its microquasar nature. The contours shown correspond to 6, 8, 10, 12, 14, 
16, 18, 20, 25, 30, 40, and 50 times 0.085 mJy beam. From Paredes et al. (2000). 



be visible to an observatory such as HESS in a matter of hours. 

The possibility of a leptonic origin for the 7-ray source cannot be ruled out in this case. The SNR 
is ranked sixth when ordered by radio flux among all entries in Green's (2000) Catalog (Table |^. With 
standard values of magnetic fields appropriate to molecular clouds (eg. Crutcher 1988, 1994, 1999), the 
radio flux that would be generated by the same electronic population producing the GeV emission would 
not overpredict the currently detected radio flux. A model such as the one presented by Bykov et al. (2000) 
could explain the EGRET source without invoking the dominance of hadronic interactions. 

10.12 7-ray source 3EG J1824-1514 - SNR G16.8-1.1 

Paredes et al. (2000) proposed that the massive star LS 5039 is part of a newly discovered microquasar 
system, and that it can be identified with the 7-ray source 3EG J1824-1514. The detection of radio jets 
and strong variability at different wavelengths (including 7-rays) supports their claim. Gamma-ray emis- 
sion from microquasars with high-mass companions has been recently discussed by Kaufman-Bernado et 
al. (2002), who show that GeV 7-rays can result from the up-scattering of UV stellar photons by the 
relativistic jet (see also Georganopoulos et al. 2002). Variability is naturally produced by the changing 
viewing angle as the jet precesses due to tidal forces from the accretion disk. 

The supernova remnant G16. 8-01.1, also within the 95 % confidence contour of the 7-ray source, has 
been recently studied by Ribo et al. (2002). Although the radio structure of the remnant is not well- 
resolved because of contamination from the partially superposed HII region RCW 164 (Rodgers et al. 
1960), but its size is ^ 30 arcmin and its total flux at 5 GHz is ~ 1 Jy (Ribo et al. 2002). The distance to 
the source is not known, but a lower limit of ~ 2 kpc has been established by Ribo et al. (2002) through 
H166a line observations of the foreground HII region. HI observations by the same authors indicate that 
the ambient density around the remnant is ~ 5 cm~^, with no evidence of interacting clouds. For a SN 
energy release E^i ~ 1 and a distance d ~ 3 kpc, we found that the expected pion-decay 7-ray flux from 
this SNR should be F{E > 100 MeV) ~ 10""^ ph cm"^ g-i (g) ^ q ^j^is g^]^^^^ ti^^^d of the 
observed flux from 3EG J1824-1514. It is possible that bremsstrahlung from SNR shell-cloud interactions 
could also generate a fraction of the observed 7-rays. 

In addition, there is one Princeton pulsar superposed on this 3EG source, but it is not energetic enough 
to be a plausible alternative to the microquasar found by Paredes et al. (see Section The apparent 
(though marginal) variability of the 7-ray source, in addition, argues weakly against a pulsar origin of the 
GeV flux. A variability analysis (Torres et al. 2001) suggests that 3EG J1824-1514 is marginally variable, 
but this is not confirmed by Tompkins's r index, which is compatible with a steady source (although the 
difference between the upper and the lower limit on Tompkins' r is large). These results are not conclusive, 
but it seems unlikely that most of the 7-ray flux could come from either the SNR or pulsar alone, especially 
when the uncertainties in distance are taken into account. 
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Figure 24: Flux history of the 3EG J1837-0423. Only for one viewing period 423.0, the source was 
undoubtedly detected, all others being only upper limits. The x-axis in the figure does not represent a 
linear scale of time. Rather, each point represent the measurement for a different single viewing period 
given in Hartman et al. (1999). 



Figure 25: The colors are CO(l-O) integrated over the velocity range 30 to 65 km/s. The solid contours 
are 4.85 GHz continuum from the survey of Condon et al. (1994). The dashed contours are the 50%, 
68%, 95%, and 99% confidence contours for 3EG J1856+0114. The white "+" marks the position with the 
highest ratio CO(2-1)/CO(1-0) as determined by Seta et al. (1998). 



The source 3EG J1824-1514 remains the best candidate for a 7-ray emitting microquasar. Future tests 
of this interesting source will surely be carried out with AGILE and GLAST. 

10.13 7-ray source 3EG J1837-0423 - SNR G27.8+0.6 

This source is variable between EGRET observations, which argues against a SNR or pulsar origin (see 
above). Indeed, 3EG J1837-0423 was detected only once, in viewing period 423.0 with 5.8 a significance 
level. In all other single viewing periods in which it was observed, only an upper limit to the flux could 
be established (see Figure . This behavior is compatible with objects presenting flares, such as AGNs. 
A microlensing model might be a plausible alternative given the absence of a strong radio emitter in the 
3EG field (Torres et al. 2002, 2003). Another possibility is a non-pulsating black hole of the sort discussed 
by B. Punsly (1998a, b and Punsly et. al. 2000). In any case, the variability clearly indicated by both the 
/ and r indices make 3EG J1837-0423 incompatible with a SNR or pulsar origin. 

10.14 7-ray source 3EG J1856+0114 - SNR G34.7-0.4 (W44) 

The possible association of the 7-ray source 3EG J1856-I-0114 and SNR G34. 7-0.4 (W44) has already been 
proposed by Esposito et al (1996) and by Dermer et al. (1997). A comprehensive review of the morpho- 
logical properties of W44 at different frequencies is given by De Jager & Mastichiadis (1996). The radio 
emission of the SNR is shell-like, but the X-ray emission is centrally peaked (Rho et al. 1994). A pulsar is 
found near the center of the 3EG source, PSR B1853-h01 (quoted as PSR J1856+0113 in Table CI). Frail 
et al. (1996) discovered a corresponding radio pulsar wind nebula, with a tail-shape pointing back to the 
center of W44. In addition, they found that the transverse velocity of the pulsar is compatible with the 
expansion speed of the radio shell. Thus, it is probable that this pulsar is the compact remnant of the 
supernova. 

Recently, two new studies concluded that the remnant is in interaction with molecular clouds: a new 
radio and optical study by Giacani et al. (1997) and a complete CO study by Seta et al. (1998). In Figure 
1251 a map of CO(l-O) integrated over the velocity range 30 to 65 km/s is shown. There are six giant 
molecular clouds, with masses between 0.3 and 3 x 10^ Mq, in the vicinity of W44. Three of them, with a 
total mass of 4.1 x IO^Mq (Seta et al. 1998) are apparently interacting with the remnant. 

The total molecular mass in the vicinity of the 3EG source in W44, specifically in the region / = 34.5 
to 34.875, b = -0.75 to -0.375, and = 30 - 65 km s'^, is 6.2 x lO'' Mq. Based on the Clemens (1985) 
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rotation curve, the main clump near the 3EG source has a kinematic distance of 2.8 kpc. The molecular 
mass is obtained based on the usual assumption that the H2 column density is proportional to the CO 
velocity- integrated intensity. The mass in the velocity-perturbed wings can be a factor ~100 less than 
the total mass of the cloud. With these values for the masses of passive targets and an CR enhancement 
factor of ~40, the entire GeV flux could, in principle, be explained by hadronic interactions, although 
an additional Bremsstrahlung component will be present also in the SNR shell-cloud interactions. The 
enhancement factor of ~40 is obtained assuming £^51 = 0.67 and a pre-shock ambient density hq = 1, but 
the conclusions are robust for reasonable variations in these parameters. 

We can compute mean densities for the well-defined clumps (lighter colors in Figure I25|) by assum- 
ing they are roughly spherical. However, the 3EG source contours enclose less intense emission that 
is spread throughout the whole W44 complex, which is about ~ 1 deg across, (49 pc at the assumed 
distance of 2.8 kpc). Thus it is reasonable to assume that the molecular gas enclosed by the 3EG 
contours is spread over about 49 pc along the line of sight as well. The mean density would then be 
n ~ 6.2 X lO^M0/[7r(9pc)2 x 49pc] = 188 H/cm'^; here, 9 pc is the 95% location contour of the 3EG 
source. This density is 6 times larger than that used by de Jager & Mastichiadis' (1997) in computing the 
hadronic emission. 

De Jager &: Mastichiadis (1997) developed a specific model for this EGRET source (at the time, 2EG 
J1857+0118) and proposed that the 7-ray radiation could be accounted for by relativistic Bremsstrahlung 
and inverse Compton scattering. One of their main motivations was to explain the hard spectral index 
of the corresponding 7-ray source, ~ —1.80, as well as the hard index of the radio source ~ —0.3, both 
difficult to reconcile with the standard Fermi first order acceleration process. The model of de Jager &: 
Mastichiadis (1997) propose the pulsar PSR B1853+01 as the source of the gamma-rays; indeed, the re- 
quired efficiency to produce all the 7-ray radiation is 13%, which appears marginally plausible, given the 
uncertainties. However, the luminosity of the pulsar wind nebula (PWN) is negligible in comparison with 
the total X-ray luminosity of W44 (Harrus, Hughes, &: Helfand 1996). The bulk of the X-ray emission from 
W44 is thermal (Jones, Smith, & Angellini 1993, Rho et al. 1994). But even though the X-ray luminosity 
of the PWN is negligible, the pulsar could have injected a significant amount of electrons in an earlier 
stage. Mastichiadis' paper Good fits of the 2EG spectrum were obtained for a range of particle density, 
whereas the magnetic field was required to be ~ 10;U G and the synchrotron cutoff frequency, Uf, ~ 10^^'^ 
Hz. A reasonable field strength could then explain the 7-ray spectrum of W44 as originating in leptonic 
processes. 

W44 is then a complicated case. While the 7-ray source may be due to SNR shock-cloud interactions 
(hadronic and leptonic) , a significant part of the 7-ray flux observed - perhaps even all of it - could be due to 
a pulsar. Application of Eq. H25|) does not allow us to discard, as in the case of SNR G347.3-0.5, a leptonic 
origin for the high-energy radiation. Here, then, future satellites and telescopes will play an essential role 
in disentangling the different possibilities. W44 was not observed to be emitting at TeV energies, where 
an upper limit has been imposed by the Whipple observatory at {F{E > 250GeV) = 8.5 x 10~^^ cm~^ 
s~^, Leslard et al. 1995), implying that there is a possible break in the spectrum from GeV to TeV. This 
can be tested by forthcoming higher sensitivity TeV-telescopes. 

10.15 7-ray source 3EG J1903+0550 - SNR G39.2-0.3 

Figure!^ shows the relative positions of 3EG J1903+0550 and the SNR G39. 2-0.3. The 3EG source is also 
coincident with two Princeton pulsars (Table 3). One of them can be readily discarded as the origin of the 
7-ray emission because of the unrealistically high value required for the efficiency. The other pulsar, PSR 
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Figure 26: Left: Relative positions of 3EG J1903+0550 and the SNR G39.2-0.3. Radio contours in step 
of 2 mK, starting from 5 niK, from the 2695-MHz map obtained with the Effelsberg 100-m single dish 
telescope (Fiirst et al. 1990). Right: CO integrated over the velocity range of the far Sagittarius arm, 
48 to 70 km s~^. Contours: 4.85 GHz continuum from the survey of Condon et al. (1994). The angular 
resolution is 7 arcmin; the contour interval is 0.5 Jy/beam, starting at 0.3 Jy/beam. The dotted circle is 
the 95% radius about the central position of 3EG J1903+0550 (Hartman et al. 1999). 



J1902+0615, lacks a measurement of the period derivative and so we cannot judge the likelihood of this 
particular association. The SNR G39.2-0.3 has been searched for OH maser emission by Koralesky et al. 
(1998), but none was detected, although, of course, this does not imply a lack of possible interaction. The 
SNR is more than 8 kpc away, which seems to be -a priori- a problem for generating the requisite flux via 
shock interactions. We find, however, that the large distance may be, at least partially, compensated for 
by the large amount of molecular material in the neighborhood. 

Caswell et al. (1975) detected HI absorption all the way up to the terminal velocity towards the SNR 
G39.2-0.3, with almost continuous strong absorption between 60 km s~^ and the terminal velocity. They 
therefore concluded that the remnant was certainly beyond the tangent point, and most likely at the far 
distance corresponding to 60 km s~^, ~9.6 kpc. Such a large distance is consistent with the high foreground 
hydrogen column inferred both by Becker & Helfand (1987) based on 21 cm absorption measurements with 
the VLA, and by Harms &: Slane (1999) based on ASCA observations. A distance of 9.6 kpc would place 
the SNR in the far Sagittarius arm, where as Figure 12^1 shows, the remnant is nearly coincident with a 
massive molecular complex. The complex is (40,59) in the catalog of Dame et al. (1986), who assigned 
the far kinematic distance based on 2 associated HII regions. The mass of this complex is estimated to be 
2.1 X W^Mq. The mass within the 95% confidence radius of the 3EG source (dotted circle in Figure I26|) 
IS even higher, 3.4 x lO^M©, b ecause the radius also includes part of another molecular complex at higher 
longitude. If part of the mass contained in the molecular complex could serve as target material for the 
relativistic particles accelerated in the SNR shock, this 3EG detection could plausibly be produced by a 
combination of Bremsstrahlung and pion decay. With the mass quoted, a CR enhancement factor of less 
than 10 is all that is needed to produce the bulk of the observed 7-ray emission. However, it is clear that 
not all of the molecular mass can be illuminated by the SNR shock front. The SNR itself is less than 8 
arcmin in size, while the 3EG source is ~1 deg in size. Only 0.1% of the molecular material need to serve 
as a target for the particles accelerated in G39.2-0.3 in order to produce the 3EG source. In this case, 
however, as in the case of 3EG J1639-4702, the enhancement factor (Ea ll5() is very large (~ 1000), as a 
result of the use of the Sedov solutions with typical values for the energy of the explosion and ambient 
unshocked density. 

There is an additional nearby SNR, G40.5-0.5 (Downes et al. 1980), which is also associated with 
the same cloud complex as G39.2-0.3 and at a similar distance. Though the center of G40.5-0.5 does 
not coincide with the 3EG source herein analyzed, it is suggesstive that the location of the 3EG source 
lies between the two SNRs, and could thus have a composite origin. (SNR G40.5-0.5 does not appear in 
Tabled due to the high ellipticity of the confidence level contours of 3EG J1903+0550.) A large region, 
—10 < & < 5, 38 < Z < 43, comprising this latter SNR as well as the 3EG source was subject of a search for 
gamma-ray emission using the HEGRA system of imaging atmospheric telescopes (Aharonian et al. 2001). 
No evidence for emission from point sources was detected, and upper limits imposed were typically below 
0.1 Crabs for the flux above 1 TeV. 
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Figure 27: CTB 87 region: The dashed contours are the usual EGRET 50%, 68%, 95%, and 99% confidence 
levels. The solid contours are taken from the 4.85 GHz continuum survey of Condon et al. (1994). The 
color is CO integrated over the range v=-65 to -50 km/s. The positions of the blazar G74.87+1.22 and of 
the WR-star WR138 are marked. 



10.16 7-ray source 3EG J2016+3657 - SNR G74.9+1.2 (CTB 87) 

Although it has been suggested that SNR G74.9+1.2 (CTB 87) may be interacting with ambient clouds 
(Huang et al. 1983, Huang and Thaddeus 1986), the coincident 3EG J2016+3657 source has been proposed 
as a counterpart of the blazar-like radio source G74.87+1.22 (B2013+370) (Halpern et al. 2001a, Mukher- 
jee et al. 2000). B2013+370 is a compact, flat spectrum, 2 Jy radio source at 1 GHz. Its multiwavelength 
properties were compiled by Mukherjee et al. (2000), and since they resemble other blazars detected by 
EGRET, make B2013+370 an interesting possible counterpart for this 3EG source. 

Optical photometry of B2013+370 shows that it is variable, providing additional evidence of its blazar 
nature (Halpern et al. 2001a). Additionally, the same authors presented a complete set of classifications 
for the 14 brightest ROSAT X-ray sources in the error circle of the 3EG source, of which B2013+370 
remained the most likely source of the 7-rays, should these come from a point-like source. The Crab-like 
supernova remnant CTB 87 is located at more than 10 kpc (Green 2000), seemingly disfavoring its shell 
interactions as the cause of the EGRET source. There are also WR stars in the field (Romero et al. 1999), 
which might produce 7-ray emission. This possibility remains to be analyzed. INTEGRAL observations 
would help in determining if there is 7-ray emission coming from the stars. 

Figure |27| shows a CO map for the CTB 87 region. One clearly defined molecular cloud appears in the 
map. The mean velocity of the molecular cloud is -57 km/s. Assuming a flat rotation curve beyond the 
solar circle, the cloud's kinematic distance is 10.4 kpc. The total molecular mass within the 95% confidence 
radius of the 3EG source is 1.7 x W^Mq. With such a high value for the molecular mass, there is still a 
chance that the hadronic or leptonic 7-ray emission may be contributing to 3EG J2016+3657. 

As in previous cases, only a precise determination of the 7-ray source position will disentangle the 
origin of this 7-ray source. Contrary to other EGRET-SNR pairs, though, this one has the particularity 
of enclosing a good candidate for an extra-galactic origin of the radiation. 

10.17 7-ray source 3EG J2020+4017 - SNR G78.2+2.1 (7-Cygni Nebula, W66) 

The SNR G78. 2-1-2.1 lies in a very complex region of the sky, where more than forty HII regions and a 
large number of shell structures exist. Recently, Lozinskaya et al. (2000) have made an in depth analysis 
of the SNR, which included new optical observations and re-analysis of archival X-ray data. We refer the 
reader to their paper for appropriate details on the SNR structure and other features, other than those 
commented here, related with the possible association with the 3EG source. The main result of Lozinskaya 
et al. is that X-ray observations lead to a self-consistent model of a young SNR at an early stage of 
adiabatic expansion into a medium of relatively low density (t = (5 — 6) x 10^ yr, no = 0.14 — 0.3 cm~'^). 

The possible association between 3EG J2020-I-4017 (with the largest signal-to-background ratio of any 
of the sources in the third EGRET catalog that are positionally coincident with shell-type supernova 
remnants) and G78. 2-1-2.1 was previously suggested by Pollock (1985), Sturner and Dermer (1995) and 
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Figure 28: The filtered radio emission at 2.7 GHz of the SNR G78. 2+2.1 is shown in black. Radio contours, 
from the 2695- MHz map obtained with the Effelsberg 100-m single dish telescope (Fiirst et al. 1990), are 
labeled in steps of 1 K in brightness temperature, starting at 3.0 K. The superposed white levels represent 
the 99%, 95%, 68%, and 50% statistical probability that a 7-ray source lies within each contour according 
to the EGRET catalog (Hartman et al. 1999). In the background, a CO(J=1-0) map of the region, 
integrated over the range v = —20 to 20 km s^^, which includes all the emission in this direction except 
for a small amount in the range —50 to —40 km s^^ which is probably associated with the Perseus Arm is 
shown. 



Esposito et al. (1996), each using the 7-ray catalog available at the time. Brazier et al. (1996) studied this 
source (2EG 2020+4626 then) and reported the discovery of a point like X-ray source, RX J2020. 2+4026, 
which lies close to the center of the remnant. If one is to assume that this source and the 3EG detection 
are related, the ratio of the gamma to X-ray fluxes is about 6000, similar to what is detected for the radio 
quiet Geminga-pulsar. This prompted Brazier et al. to suggest that 3EG J2020+4017 can indeed be a 
new Geminga-like radio quiet pulsar, something which would be in tune with the hard spectra, the low 
variability index, and the absence of 'em' classification in the 3EG Catalog (see Table This possible 
association remains, then, very suggestive. However, Brazier et al. have reported no evidence for pulsations 
in the range 13.3 > /(Hz)> 2, F < 10^^^ s~^. Also, no radio pulsar is known to be superposed to the 3EG 
contours (see Table Ej). The absence of pulsar 7-ray radiation can well be due to a statistical limitation of 
the EGRET data. Analysis of GLAST (LAT) performance shows that periodicities should be detectable 
if present in any of the known low-latitude EGRET sources (Carramihana 2001, Thompson et al. 2001). 

In particular. Brazier et al. considered that the absence of em (i.e. extended) classification in the 3EG 
Catalog was definitive in disregarding the SNR as a possible site of cosmic ray acceleration: the SNR is 
about 1 deg in size, and it should be visible as extended by EGRET if the ambient matter were uniformly 
distributed. This, however, may not be the case if the shock accelerated particles interact and emit 7-rays 
most intensely at a localized nearby concentration of molecular material. 

Yamamoto et al. (1999) have studied the possible SNR-cloud interaction here and observed a very 
high CO(J=2-1)/CO(J=1-0) intensity ratio (~ 1.5), very suggestive of an interacting cloud, is observed at 
(/, 6)=(78,2.3). Interestingly, this position exactly coincides with that of the 3EG source (see Figure 05)) . 
After re-analyzing the same set of data used by Yamamoto et al. (1999), we agree with their result, but 
have also found another high ratio CO(J=2-1)/CO(J=1-0) at an adjacent position (/ = 77.875, b = 2.25), 
and yet another moderately high ratio (~ 1) at another adjacent position: {l,b) = (78.00,2.25). These 
high ratios coincide nicely with the 7-ray source and with a fairly well-defined CO cloud, and make a good 
case for the interaction between the SNR and the cloud. 

Assuming a distance of 1.7 kpc (Lozinskaya et al. 2000), the molecular mass of this cloud is 4700 Mq 
(calculated over the range / = 77.75 - 78.125 deg, b = 2.00 - 2.375 deg, and v = -10 to km s"^). This 
mass is rather uncertain since W66 lies in the so-called Cygnus X region of the Galaxy where the CO 
emission is very complex and strong, probably because we are viewing the Local spiral arm tangentially. 
It is hard to judge what amount of CO might be associated with the SNR and what amount is seen in 
projection. The Local arm emission is mainly in the range —20 to 20 km s~^, and the Perseus Arm is seen 
at more negative velocities. Kinematic distances, particularly at low velocity, are very unreliable because 
our line of sight is almost tangent to the solar circle. 
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Using the data obtained for the mass of the cloud in the vicinity of the 3EG source and apparently 
in interaction with the SNR, as well as the mean value for the unshocked density (no = 0.22 cm~^) and 
distance (from Lozinskaya et al. 2000), it is possible to explain the 7-rays by hadronic processes. The use 
of Eq. (|25j) with plausible values of the magnetic fields in molecular clouds (Crutcher 1988, 1994, 1999) 
rules out leptonic processes as the source of most of the radiation. 

We cannot discard, however, a composite origin for the 7-rays, with part of the radiation coming from 
the putative 7-ray pulsar proposed by Brazier et al. (1996) and part from cloud interactions with the 
nucleonic component of freshly accelerated CRs. 

10.18 An example beyond |6| > 10: 7-ray source 3EG JOOlO+7309 - SNR 
G119.5+10.2 (CTA 1) 

The SNR-EGRET connection that we have explored so far was based on the sample shown in Table 1, 
constructed using 3EG sources within 10° of the Galactic plane. Although most SNRs fall into this latitude 
range, there is the possibility of finding a few nearby SNRs related with 7-ray sources at higher latitudes. 
An example of such is briefly described in this section. 

CTA 1 is a shell-type SNR, but its shell is incomplete and broken-out towards the NW. This breakout 
phenomenon may be caused by more rapid expansion of the blast wave shock into a lower density region 
toward the NW. HI observations supported this interpretation (Pineault et al. 1993, 1997). Since CTA 
1 is located at a relatively high latitude (10 deg) and is nearby (1.4ib0.3 kpc, Pineault et al. 1993), it 
has a large angular size (90 arcmin), little foreground or background confusion, and it can be observed 
at exceptionally high linear resolution. The age of the SNR was estimated to be 10^ yr by Pineault et 
al. (1993), but it could be younger by a factor of 2 (Slane et al. 1997). CTA 1 was subject of intense 
observational campaigns in the past years. There have been both ROSAT and ASCA X-ray observations 
(Seward et al. 1995, Slane et al. 1997), as well as optical, infrared and radio (see Pineault et al. 1997 and 
Brazier et al. 1998 for a review). 

The ROSAT observation confirmed that CTA 1 belongs to the class of composite SNRs, which show a 
shell-type morphology in the radio band and are centre-filled in X-rays. Five point sources were detected 
with ROSAT, one of which was found to coincide with the EGRET source (at the time of the analysis, 2EG 
J0008+7307). ASCA data later revealed that this source, named RX J0007. 0-1-7302, has a non-thermal 
spectrum, suggesting that it is the pulsar left from the supernova explosion (Slane et al. 1997). Optical 
observations were carried out by Brazier et al. (1998), with a 2.12-m telescope, but no object was found 
within the positional error box of the X-ray source. This allowed an upper limit to be set on the optical 
magnitude of any counterpart to the putative pulsar. 

The 3EG JOOlO-l-7309 is a non-variable source under the / and r schemes, and has a hard spectral 
index of 1.85ib0.10, compatible with those of the Vela pulsar. This source was also detected in the Second 
EGRET Catalog, but with a shifted position which made it coincide with a nearby AGN - which was at 
the time proposed as a possible counterpart (Nolan et al. 1996). This AGN is not coincident with the 
3EG source and can not be considered a plausible counterpart any longer. Based on positional coincidence, 
on the hard spectral index, and on physical similarities between the Vela pulsar and RX J0007. 0-1-7302, 
Brazier et al. (1998) proposed that the 3EG source and this X-ray source were related. For an assumed 1 
sr beaming, the observed 100-2000 MeV flux corresponds to a luminosity of 4 x 10^^ erg s~^, compatible 
with other 7-ray pulsar detections (see Table |H| • 
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Figure 29: The Capricornus SNR uncovered by the hkely associated high-energy emission. Left: Back- 
ground filtered radio emission at 408 MHz (radio data from Haslam et al. 1982) of the region surrounding 
three 3EG sources, whose contours are marked. Middle: Radio map at 2.3 GHz (data from Jonas et al. 
1998). The use of both maps shows that the SNR is a non-thermal radio source. Right: A map of the 
integrated column density of HI (velocity interval -3 to -|-5 km s~^). Label units are 10^^ atoms cm~^. Use 
of this map shows that small enhancement factors could produce all three EGRET sources. From Combi 
et al. (2001). 



Although a thorough CO investigation has not yet been done for this SNR, current data from HI 
observations do not support the presence of very dense and massive molecular clouds in its neighborhood 
(Pineauh et al. 1997). It seems CTA 1 and 3EG J0010-h7309 might be related only through the compact 
object left by the latter (RX J0007. 0-1-7302). A better localization of the 7-ray source with AGILE and 
GLAST will certainly test this suggestion. 

11 SNRs discovered by their likely associated high-energy radiation 

The SNR catalog compiled by Green (2000) is by no means complete. A large number of low surface 
brightness remnants remain to be discovered, hidden in the diffuse non-thermal continuum emission pro- 
duced by the diffuse component of the cosmic-ray electrons in the Galaxy. In recent years, the application 
of filtering techniques in radio continuum data has lead to the detection of several new SNR (e.g. Duncan 
et al. 1995, 1997, Jonas 1999, Combi et al. 1999, Klothes et al. 2001). The use of unidentified non- variable 
7-ray sources as tracers of interacting remnants can result in new SNRs being discovered by their (plausibly 
associated) high-energy emission. 

As it was emphasized in the previous sections, not all SNR generate observable 7-rays (say, at EGRET 
sensitivity) at typical Galactic distances. A second ingredient is necessary: a target, i.e., a dense medium 
such as a molecular cloud. If a SN explodes in a cloudy medium then more than a single cloud could be 
illuminated hy p — p collisions and thus multiple 7-ray sources can emerge. Clusters of steady unidentified 
7-ray sources could trace these situations and lead to the discovery of new, very extended SNRs of low 
radio surface brightness. 

This approach of considering EGRET sources as tracers of SNRs has been applied by Combi et al. 
(1998, 2001), leading to the discovery of two new SNRs. The basic technique consists of making HI line 
observations in the direction of clearly non-variable 7-ray sources. If some well-defined but small clouds 
(M ~ 10^ — 10"^ -^0) are found within the 95% contours of the 3EG sources at velocities that correspond, 
according to the galactic rotation curve, to distances of less than 1 kpc, then large-scale (i.e. several de- 
grees) radio continuum observations are carried out. These observations aim at detecting, through image 
filtering techniques, large SNRs of low surface brightness. Observations at more than a single frequency are 
necessary in order to determine if the radio spectral index is non-thermal, as expected from such remnants. 

The SNR candidates discovered by this procedure are G327-12.0, located in the ARA region (Combi et 
al 1998) and G06. 5-12.0, located in Capricornus (Combi et al 2001): see Figure EHl The first one appears 
to be responsible for the 7-ray source 3EG J1659-6251, whereas the second could be related to a cluster of 
three 3EG sources: 3EG J1834-2803, 3EG J1847-3219, and 3EG1850-2652. HI clouds have been found at 
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the positions of all these sources except for 3EG J1847-3219. Both remnants are thought to be so close that 
cosmic-ray enhancement factors in the range 5 < feg < 45 are sufficient to explain the 7-sources through 
vr'^-decays alone. Their low radio surface-brightnesses rule out electronic Bremsstrahlung as the source of 
the GeV flux. 

We expect that the use of filtering techniques in interferometric radio observations could lead to the 
discovery of more such remnants in the near future (see, eg., Klothes et al. 2001), and with the improved 
capabilities of the next generation satellites, 7-ray emission from far more distant interacting SNRs could 
be detected. 

12 SNRs and their neighborhoods as TeV sources 

As can be seen in Table |71 several observations of SNRs such as W44, W51, 7-Cygni, W63 and Tycho's 
SNR, selected because of their possible association with molecular clouds and/or EGRET sources, have 
only produced TeV emission upper limits (Buckley et al. 1998). This, as discussed below, could indicate 
spectral cutoffs or breaks in the GeV- TeV energy range. For instance, the required differential source spec- 
trum would have to steepen to ~ E~^-^ for 7-Cygni in order to escape detection at TeV energies (Fegan 
2001). 

Very recently (Aharonian et al. 2002a), the HEGRA system of imaging atmospheric Cherenkov tele- 
scopes reported a survey of one quarter of the Galactic plane (—2 < I < 85). TeV gamma-ray emission 
from point sources and moderately extended sources (diameter less than 0.8 deg), including 86 known pul- 
sars (PSR), 63 known supernova remnants (SNR) and 9 GeV sources, were searched with negative results. 
Upper limits range from 0.15 Crab units up to several Crab units, depending on the observation time and 
zenith angles covered: no TeV source was detected above 4.5cr in a total observation time of 115 h. At 
the same time, a search for point sources of radiation above 15 TeV has been conducted with HEGRA 
AIROBICC array (Aharonian et al. 2002b), but only flux upper limits of around 1.3 times the flux of Crab 
nebula were obtained for candidate sources (including SNRs), depending, again, on the observation time 
and zenith angles covered. 

Positive detections, however, already exist. SN1006 (Tanimori et al. 1998) and RXJ1713. 7-3946 (Mu- 
raishi et al. 2000, Enomoto et al. 2002) were detected by the CANGAROO telescopes. Observations of 
SN1006 in 1996 and 1997 show a significant excess from the NW rim of the SNR. The excess is consistent 
with the location of non-thermal X-rays detected by ASCA (Koyama et al. 1995). Theoretical modeling 
for SN 1006 has also been performed (e.g. Berezhko et al. 2003). 

Cassiopeia A was recently announced as TeV source by the HEGRA collaboration (Aharonian et al. 
2001c). 232 hours of observations yield an excess at the 4.9cj level, and a flux of F — 5.8 i 1.2gtat i 
2syst X 10~^^cm~'^s~^ at (E> 1 TeV). The origin of the 7-rays has been recently discussed by Berezhko 
et al. (2003). Cas A has already been associated with a bright source of hard X-rays which indicates a 
population of non-thermal electrons with energies up to 100 TeV (Allen et al. 1999b). Nevertheless, Cas 
A is not a 3EG source, but this could be due only to the small EGRET sensitivity. 

There have also been observations of unidentified 7-ray sources whose identifications are, perhaps, more 
tentative. These observations were made with the Whipple 10-m telescope (Buckley et al. 1997) Two of 
them were related with the 3EG sources in our sample, and details are provided in Table [Tj 
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Table 7: TeV Observations of SNRs. Plerions are shown in the first panel and shell-type SNRs in the 
second. The third panel shows results for two 3EG sources coincident with SNRs for which there have also 
been TeV observations. Partially adapted from Fegan (2001) and Mori (2001). 







TTli TV / T Tnnpv T .1 TTi 1 ■f r\v aT\cir''f viitvi 

lU-rt./ kJJJJJcl iJililit KJL OjJcL/tl Ulli 

in — 11^™— 2„ — 1 ^„ in— 11 p™— 2 „ — 1 Tp^aV"^ 
lU Cflt b Ul lU Cill o V 


IN dilllt. 




A 11 Tf^A/ (^Vici<=»i^vF»^"riTi<^c: 






\jLaiU INCLILiid 




7 4oonpv^ 

( .W \^ rt\J\J\JK^ V 1 


CANOAROO 






PSR 17(16-44 


60 


15 ('>lTf'V') 

\J • X tj V ^ X X.C V 1 


V CXd X LlXodX 


116 


n 26 (Fj/2 TpV")"^-* TpV~^ 


' ' U.X J.XCI.XXJ. 






PSP! 1706-44 


10 


1 2 f>300C;oV') 




H.7~) 




D ANH AROO 






T?X T 1 71 7 "^QAf^ 
xl-A-J lilo.r-oy40 


00 


u.oo \^i--o lev J 


RX7 1718 7-3946-^ 


32 


53 1 63 ± 15 ± 321 R-2-84±o.i5±o.20 




"^4 


U.4:U [^^•t -Lev ( 


VV ZrO 


oo 


*\U.OO \^ O -Lev I 


HEGRA 






Cas A 


232 


0.058 (> 1 TeV)^ 


7-Cygni 


47 


<1.1 (>500GeV)'= 


Monoceros 


120 




Durham 






SN1006 


41 


<1.7 (>300GeV) 


Whipple 






Monoceros 


13.1 


<4.8 (>500GeV) 


Cas A 


6.9 


<0.66 (>500GeV) 


W44 


6 


<3.0 (>300GeV) 


W51 


7.8 


<3.6 (>300GeV) 


7-Cygni 


9.3 


<2.2 (>300GeV) 


W63 


2.3 


<6.4 (>300GeV) 


Tycho 


14.5 


<0.8 (>300GeV) 


CAT 






CasA 


24.4 


<0.74 (>400GeV) 


J2016-f-3657 


287 


5.8^* 


J2020-F4017 


513 


0.990^^ 



"A different definition of Energy Threshold is used. ''Evidence for emission at the 4.9(t level (Piihlhofer et al. 2001). 
^Limits converted from Crab units using flux of Hillas ct al. (1998). ''integral Flux Above 400 GeV. «Not yet 
reported. ^ Observations made with CANGAROO-II (Enomoto et al. 2002). 
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The Crab Nebula has been detected by several Cherenkov observatories and it is often used to provide 
a check on the calibration of new instruments (e.g. Aharonian ct al. 2000). STACEE and CELESTE, 
on which we comment below, have published significant detections of the Crab Nebula in an energy range 
lower than that obtained by others ground-based instruments, E > 190 ± 60 GeV (Oser et al. 2000) and 
E > 50GeV (De Naurois et al. 2001), respectively. The Crab's energy spectrum between SOOGeV and 
50TeV has been well established (sec below), and it steepens with energy (Hillas et al. 1998, Aharonian 
et al. 2000). No pulsation has yet been seen (Gillanders et al. 1997, Burdett et al. 1999, Aharonian et al. 
1999). 

The CANGAROO team has detected the pulsar PSR 1706-44 in 60 hours of observations (Kifune et 
al. 1995). They have also detected the Vela pulsar at the 6a level, based on ~ 120 hours of observation. 
Observations by the Durham group (Chadwick et al. 1997) also confirmed these detections. In the case of 
Vela, the VHE signal, which is offset from the location of the pulsar by 0.14°, is thought to originate from a 
synchrotron nebula, powered by a population of relativistic electrons which were created in the supernova 
explosion and which have survived since then due to the low magnetic field in the nebula. 

13 Concluding remarks 

The coming years will be exciting times again for 7-ray astronomy - after the unfortunate forced demise 
of the Compton Observatory. INTEGRAL, AGILE, MAGIC, and the new stereo-IACTs are or should 
be on-line soon; and in the longer term, GLAST will surely bring about another revolution, answering 
existing questions, such as the one reviewed here, and posing further challenges. Ultimately, it will take 
a sensitive and high spatial-resolution MeV-GeV detector such as GLAST working in close consort with 
the ground-based TeV and radio telescopes to address the Galactic cosmic ray origin problem: the single 
most definitive test that SNRs (or any other putative sources) are accelerating CR nuclei would be the 
statistically significant detection of the signature neutral pion 7-ray hump centered at 67.5 MeV (in log£'-y), 
as has been seen already by EGRET in the diffuse 7-ray background (Hunter et al. 1997). Though such 
a detection, by itself, only proves the existence of lower energy (~1 GeV/n) nuclei, the detection allows 
to normalize the hadronic vs. electronic contributions in a model-independent fashion at those 'lower' 
energies. Thus, such a detection together with the extension of the spectrum into the TeV regime, and a 
multiwavelength spectrum inconsistent with an electronic origin, ought to be sufficient evidence to conclude 
that nuclei are being accelerated by SNRs. It remains to be seen if this will cleanly be borne out by the 
data. 

In any case, the SNR-EGRET connection looks stronger than ever. Several cases for a physical associa- 
tion, based on the analysis of gamma-ray data and the molecular environment of the SNR, has been shown 
promising and a definite conclusion about the origin of cosmic rays in SNR shocks seems to be reachable 
in a nearby future. Very recently, Erlykin and Wolfendale (2003) have shown, in addition, that for some 
nearby SNRs, and even when there is no gamma-ray signature, the idea of cosmic ray production in the 
shocks cannot be discarded since there could be a possible evacuation of ambient gas by the stellar wind 
of the progenitor star, or by the explosion of a nearby earlier supernova. The case by case analysis shows, 
moreover, that it is at least plausible that EGRET has detected distant (more than 6 kpc) SNRs. There 
are 5 coinciding pairs of 3EG sources and SNRs for which the latter apparently lie at such high values of 
distance (disregarding those related with SNRs spatially close to the galactic center). For all these cases, 
we have uncovered the existence of nearby, large, in some cases giant, molecular clouds that could enhance 
the GeV signal through pion decay. It is possible that the physical relationship between the 3EG source 
and the coinciding SNR could provide for these pairs a substantial part of the GeV emission observed. 
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Table 8: Instrumental parameters of some of the forthcoming satellite missions in the MeV-GeV range. 
EGRET data (first panel) is shown for comparison. The second panel corresponds to AGILE-GRID and 
the third to GLAST-LAT. 
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This does not preclude, however, composite origins for the total amount of the radiation detected. Some 
of these cases present other plausible scenarios. AGILE observations, in advance of GLAST, would greatly 
elucidate the origin for these 3EG sources, since even a factor of 2 improvement in spatial resolution would 
be enough to reject the SNR connection. 

Kilometer-scale neutrino telescopes have also proposed as viable detectors of hadronic cosmic ray sources 
(Halzen & Hooper 2002, Anchordoqui et al. 2003), and will be a welcomed addition to the arsenal of space- 
and ground-based detectors that ought to be lined up by the time the large-scale neutrino telescopes are 
functional. Surely, with all this instrumentation focused on the problem we will finally be able to test 
Shklovskii's suspicion that "it is possible that ionized interstellar atoms are accelerated in the moving 
magnetic fields connected with an expanding [SNR] nebula." (Shklovskii, 1953). 

14 Appendix: Reviewing the prospects for the forthcoming GeV satel- 
Htes 

14.1 INTEGRAL 

The International Gamma-ray Astrophysics Laboratory (INTEGRAL) is now in orbit. It has main sci- 
entific instruments, named SPI (Spectrometer of INTEGRAL) and IBIS (Imager onBoard the Integral 
Satellite). SPI will work in the range 20 KeV - 8 MeV and will perform high resolution spectroscopy; it 
can measure 7-ray line profiles with an accuracy of 2 keV. The sensitivity of IBIS lies in the range 10 keV 
- 10 MeV, and it will focus on achieving good angular resolution (12 arcmin). INTEGRAL will also have 
three monitors, the twin JEM-X's (3 keV - 35 keV) and a CMC in the optical band (500-600 nm). The 
JEM-X's will have an angular resolution of 3 arcmin and a 4.8° fields of view. The CMC will have a pixel 
resolution of 16.6 arcsec and a 5° x 5° field of view (Schonfelder 2001). 

INTEGRAL could identify compact objects with good spatial resolution that are likely counterparts of 
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EGRET sources, and will help confirm or reject interpretations based, for instance, on microquasars and 
7-ray blazars. The INTEGRAL galactic plane exposure will help to clarify the confused region near Vela. 
Another 1 Ms exposure, part of the satellite Open Program, will focus on the Garina region. That part 
of the sky is populated with five EGRET sources, two of which were analyzed above in the case-by-case 
study: 3EG J1102-6103 and 3EG J1013-5915. As an example of its potential, we will briefly discuss what 
INTEGRAL can do in establishing the nature of these sources. 

The former case, 3EG Jl 102-6103, was briefly mentioned in the corresponding section above. Although 
a hadronic origin for the 7-ray emission is possible for this source, it could also be the result of inverse 
Gompton up-scattering of UV photons by electrons accelerated in the winds of one or several Wolf-Rayet 
stars (particularly in the region where the winds of W39 and WR38B collide) . If the latter explanation is 
correct, INTEGRAL should see a source where the winds collide. 

In the latter case, the source 3EG J1013-5915 is likely mostly produced by PSR J1013-5915 (Gamilo 
et al. 2001). The JEM-X monitor, in particular, could then probe the putative non-thermal emission 
from this region, especially the tail of the synchrotron spectrum and, in this way, help to determine the 
high-energy cutoff of the electron population in the source. 

In addition, there have been suggestions that INTEGRAL could detect the 0.511 MeV line from giant 
molecular clouds (see e.g. Guessoum et al. 2001). Giant molecular clouds are typically surrounded by HII 
regions, ionized in an uncertain fraction. If cosmic-rays are diffusing into the cloud, the cores, too, may be 
ionized. If a sufficient density of cosmic-rays is present, they could excite GNO nuclei as well as lead to the 
production of nuclear gamma-rays lines and positron-production, the latter being emitted by radioactive 
nuclei. The fluxes and line separations predicted are on the verge of detectability by INTEGRAL (Gues- 
soum et al. 2001), but it will be an interesting arena to explore with forthcoming observations, particularly 
for those giant molecular clouds closer to Earth (regrettably not the ones we find superposed with the SNRs 
in Table 1). 

Thus, although it is expected that no direct observations with INTEGRAL can prove that cosmic- 
rays are accelerated in SNR shocks (the energy range of the relevant phenomena, around 100 MeV and 
beyond, is out of the INTEGRAL energy band), they can be used to explore alternative explanations 
for the unidentified EGRET sources, and so are important to get an overall picture of the SNR-EGRET 
source connection. Nuclear gamma-ray line signatures in the 0.1-10 MeV band will certainly also provide 
important hints regarding the GR acceleration processes in the Galaxy. 

14.2 AGILE 

AGILE (Astro-rivelatore Gamma a Immagini Leggero), is expected to be launched in 2003 (for a recent 
review see Tavani et al. 2001). AGILE will have a very large field of view, covering at one time approxi- 
mately 1/5 of the sky at energies between 30 MeV and 50 GeV. The angular resolution will be a factor of 
two better than that of EGRET, specifically due to the GRID instrument, whose parameters are given in 
Table IHl However, the sensitivity for point-like sources will remain comparable to that of EGRET. AGILE 
will also have detection and imaging capabilities in the hard X-ray range provided by the Super-AGILE 
detector. The main goals of Super-AGILE are the simultaneous 7-ray and hard X-ray detection of astro- 
physical sources (which was never achieved by previous 7-ray instruments), improved source positioning 
(1-3 arcmins, depending on intensity), fast burst alert, and on-board triggering capability. 

AGILE is thus very well suited for studying compact objects, particularily those presenting 7-ray 
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Figure 30: Models of the neutral pion decay, non-thermal Bremsstrahlung (NB), inverse Compton (IC), 
and pulsar (PSR) 7-ray spectra associated with W66. The sum of the three cosmic-ray components (vr" -|- 
NB -|- IC) of the shell and the sum of all four components (vr^ -|- NB -|- IC + pulsar emission) are shown. 
EGRET spectral data is also included. From Allen et al. (1999). 

variability or pulsed emission. AGILE will search for pulsed 7-ray emission from all recently discovered 
Parkes pulsars coincident with EGRET sources (D'amico et al. 2001, Torres et al. 2001d, Camilo et al. 
2001), so establishing their contribution to the EGRET 7-ray flux. Unpulsed 7-ray emission from plerionic 
SNRs, and search for time variability in pulsar- wind nebula interactions will also be possible targets for 
AGILE. Finally, AGILE will be essential in assessing the possible existence of new populations of variable 
7-ray sources in the galaxy, such as non-pulsating black holes, X-ray binaries, and microquasars. In the 
case of 3EG J0542-I-2610, for instance, AGILE could test the hypothesis that the 7-ray emission is produced 
in A0535-I-26 (Romero et al. 2001). The key prediction of this model is anti-correlation of the X-ray and 
7-ray emissions. 

14.3 GLAST 

The Large Area Telescope (LAT) on the upcoming Gamma-ray Large Area Space Telescope (GLAST) 
mission will be suitable for studying the relationship between 7-ray sources and SNRs. The instrument's 
parameters are given in Table IHl (Michelson 2001). GLAST, expected to be launched in 2006, will explore 
the energy range from 30 MeV to greater than 100 GeV with 10% energy resolution between 100 MeV and 
10 GeV. GLAST uniquely combines high angular resolution with superb sensitivity, and has a moderate 
effective area. Sources below the EGRET threshold (~ 6 x 10^® photons cm~^ s~^) will be localized to 
arcmin scales. This clearly will improve our understanding of the SNR shock acceleration of cosmic-rays 
and hadronic 7-ray production. 

In what follows, we provide a brief example of GLAST capabilities as applied to the SNR 7-cygni 
(W66), courtesy of Seth Digel and NASA (see also, Ormes et al. 2000 astro-ph/0003270 from which this 
example is adapted). Allen et al. (1999) studied what information could be obtained from GLAST obser- 
vations based on a 1-year all-sky survey assuming that 60% of the 7-ray flux was produced by a pulsar at 
the location proposed by Brazier et al. (1996). This pulsar was assumed to have a differential photon index 
of 2.08, that of the coincident source, 3EG J2020-I-4017. The remainder of the photon flux was assumed to 
come from relativistically accelerated particles through both leptonic and hadronic processes - the large 
dominance of the latter in agreement with what we found in the corresponding section above. The position 
of the molecular cloud assumed in the simulations corresponds well with the CO clouds actually found - 
see Figure 0H|) . 

Allen et al. additionally assumed that the electron and proton spectra of W66 have the shape specified 
by Bell (1978, see their Eq. 5), with a common relativistic spectral index of F = 2.08. The normalization 
of the electron spectrum was determined from the radio data, by assuming that the magnetic field strength 
is 100 fiG. The normalization of the proton spectrum, instead, was determined by assuming that the total 
number of non-thermal electrons is a factor of 1.2 larger than the number of non-thermal protons. 

The spectral results of the simulations are shown in Figure 1301 It shows the 7-ray spectra produced by 
the putative pulsar, by the decay of neutral pions, as well as by the leptonic processes: Bremsstrahlung 
radiation of the electrons, and inverse Compton scattering of electrons on the cosmic microwave background 
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Figure 31: Comparison between the observed EGRET data and a GLAST simulation for the W66 SNR 
region. The large circle shows the position and extent of the radio shell of the SNR. GLAST will be able to 
localize distinctively the position of the putative 7-ray pulsar (marked here as the X-ray source) proposed 
by Brazier et al. as well as the SNR shock interacting with the molecular cloud, should both contribute 
as assumed to the 7-ray flux observed. Courtesy of S. Digel and NASA; adapted from Ormes et al. 2000 
astro-ph/0003270 



radiation. The latter three of these four spectra were obtained using the 7-ray emissivity results of Gaisser, 
Protheroe, & Stanev (1998), and Baring et al. (1999). The non-thermal Bremsstrahlung and neutral pion 
spectra were obtained assuming that the average density of the material with which the cosmic-rays interact 
was no = 190 atoms cm~'^, which should be compared with the actual density we found, of 188 atoms 
cm~^ (see Section [10.171 Perhaps most illustrative is Figure I5TI which demonstrates that the resolution of 
GLAST will allow one to distinguish the contribution of the pulsar from that of the interacting molecular 
cloud. Should a picture like Figure lSTl be the result of an actual GLAST observation the proposed composite 
origin for 3EG J2020-F4017 would be proved. 

15 Appendix: Future TeV telescopes and their look at SNRs - Adapted 
from Petry (2001) 

At higher and higher energies the sky appears darker and darker. There are just 57 sources detected be- 
tween 1 and 10 GeV (Lamb & Macomb 1997) and higher frequencies remain largely unobserved. In order 
to reach the highest energy 7-rays, new ground-based telescopes are being built and older ones upgraded. 
A very brief description of them, and of their impact on the SNR-7-ray source association problem is given 
here. 

One approach to constructing high-energy 7-ray detectors is to use existing solar farms, which have 
fields of large heliostats focusing sunlight on a central tower; such facilities lie unused at night. The arrival 
of the Cerenkov wavefront at groups of heliostats is precisely measured, and this information is used to 
differentiate 7-rays from cosmic-ray primaries (Ong 1998). STACEE (Chantell et al. 1998), CELESTE 
(Smith et al. 1997), Solar-2 (Tiimer et al. 1999), and GRAAL (Arqueros et al. 1999) are ah examples of 
such facilities. 

MILAGRO (Sinnis et al. 1995) and Tibet HD (Amenomori et al. 1999) are examples of air showers 
detectors. These arrays can in principle operate 24 hours a day and are expected to achieve a larger energy 
range, comparable to the next generation imagers discussed below. 

Imaging Atmospheric Cherenkov Telescopes (lACTs) are ground-based 7-ray detectors using the at- 
mosphere as a tracker and calorimeter and having good all-around performance, both in sensitivity and 
angular resolution, for sources above 100 GeV. They typically have point spread functions with 9qs < 0.16°. 
The number of detected primary gamma-photons will typically be > 100, and therefore source locations, 
computed as 9qs/^/N, can reach arcmin accuracies. Generically, the possibility of separating two nearby 
point sources can be realized if their angular distance is > 36*68 • Some forthcoming lACT telescopes in the 
TeV regime are detailed in Table El 

The successful introduction of stereo imaging technology in the HEGRA instrument has paved the way 
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Table 9: Next-generation Imaging Cherenkov Telescopes and their main characteristics. The locations of 
the sites chosen for these observatories and their estimated minimum energy thresholds are given. Roman 
numbers behind some of the project names denote the project phases. The Type refers to the number of 
individual telescopes and diameter of their mirror dish. The last column is the predicted energy threshold 
for 7-ray photons at zenith angle 1? = 0°. (From Petry 2001). 
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for bigger and more complex instruments. The unprecedented accuracy in air shower stereo reconstruction 
of the HESS telescopes in Namibia will allow one to obtain an angular resolution of ~ 0.1° with a source 
location accuracy of only 10 arcseconds (Konopelko 2001). After 50 hours of observation of a point source, 
HESS will detect a flux of ~ 10~^^ ph cm~^ s~^ at £^ > 100 GeV. In the case of extended sources of 
angular size 9, the flux detected will be ~ (0/0.1 deg) x 10"^^ ph cm~^ s~^ in the same energy range. It 
is clear that such a powerful instrument will provide answers to many of the pending questions mentioned 
in this review, at least for the case of southern SNRs. 

A natural next step in the development of stereo imaging arrays for TeV astronomy is to place one of 
these systems at high altitude, where it could detect much lower-energy gamma-rays. It has been estimated 
that at an altitude of 5 km, a threshold of 5 GeV might be achieved (Aharonian et al. 2001b). These 
systems, then, might replace orbital observatories in the GeV band in the foreseeable future. 

In an attempt to establish an agenda for the new TeV telescopes, Petry (2001) and Petry Sz Reimer 
(2001) have analyzed the observational possibilities for each of the unidentified and tentatively identified 
EGRET sources. Several problems have to be taken into account. The first is that of the apparent spectral 
steepening of several 3EG sources between the GeV and the TeV band (see Reimer and Bertsch 2001). 
This could be why strong sources such as 7-Cygni, IC443, W28, and CTA 1 were not observed in the TeV 
band (eg. Buckley et al. 1998). In order to get an estimate of the impact of this effect on any 3EG source, 
Petry (2001) proposed to use the same spectral steepening as that of the Crab Nebula. The differential 
spectral index Fo.iGeV of the Crab Nebula at 0.1 GeV is Fo.i = 2.19 it 0.02 (Hartman et al. 1999) while it 
was measured by the Whipple telescope at 500 GeV to be FsooGeV = 2.49 ± 0.06 ± 0.04 (Hillas et al. 1998). 
The latter authors showed that the steepening towards higher energies can be described by an increase in 
spectral index of 0.15 per decade of energy. Petry (2001) showed that this would then imply, if a similar 
situation is the case for the unidentified EGRET sources, that 

F{E > x[Ge\]) = Fo ■ 10"" • 10-("+0-i5) . ( x ^-(a+o.30) 
= Fo • 10-("-°-^5) • 
for 10 < X < 100 

(39) 

F{E > x[GeY]) = Fq ■ 10"" • lO-^^+O-^^) . io-("+o.30) . (_x_^-(a+o.45) 

= Fo • 10-("-0-'^5) . ^-(a+0.45) 

for 100 < X < 1000 

where Fo is the integrated flux above 0.1 GeV and a is the flux spectral index (obtained by subtracting 



54 



1.0 from the differential spectral index T, quoted for instance in Table 1) of a given source taken from the 
Third EGRET Catalog. An extrapolation without taking into account possible breaks in the spectrum 
could give a large over-estimation of the high-energy flux. Although not valid for all sources, the previous 
scheme can be used to give an idea of the expected flux from interesting 3EG detections (Petry 2001, Petry 
& Reimer 2001). Nonetheless, not only the spectral cutoffs, but other technical problems affect actual 
observation. ^ 

Taking these considerations into account it is possible to determine which of the 3EG sources analyzed 
here might be detected by the next generation lACTs. 12 out of the 19 SNR-EGRET cases listed in Tabled 
can be considered likely candidates, from a practical point of view. The complete results for all 3EG sources 
superposed to SNRs that might be detected in less than 50 hours is compiled in Table El Columns are 
as follows: 3EG name, radius of the 95% confidence contour, spectral index at 0.1 GeV, all from Hartman 
et al. (1999), minimum zenith angle, lACTs minimum energy threshold, expected flux at the minimum 
energy threshold, the integral spectral index at the lACTs minimum energy threshold for this source, the 
observation time to obtain a detection with 5 a significance (in brackets are the corresponding data if the 
spectrum at 0.1 GeV was steeper by one standard deviation). A ^-mark indicates that the detection is 
photon flux limited: the observation time was increased such that 100 photons are detected. 

Many of the SNR-EGRET pairs studied will be primary candidates for one or several of the forthcoming 
lACTs. A combination of TeV and GeV observations, together with an understanding of the molecular 
material distribution of the region, will be crucial in determining the nature of several of the unidentified 
EGRET sources analyzed here. 
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Table 10: TeV observability of EGRET-SNR pairs, adapted from Petry (2001). 
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